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We present details of simulations for the light hadron spectrum in quenched QCD carried out 
on the CP-PACS parallel computer. Simulations are made with the Wilson quark action and the 
plaquette gauge action on lattices of size 32 3 x 56 — 64 3 x 112 at four values of lattice spacings in the 
range a m 0.1-0.05 fm and the spatial extent L s a w 3 fm. Hadronic observables are calculated at five 
quark masses corresponding to mps/mv « 0.75-0.4, assuming the u and d quarks being degenerate 
but treating the s quark separately. We find that the presence of quenched chiral singularities is 
supported from an analysis of the pseudoscalar meson data. Physical values of hadron masses are 
determined using m^, m p and ttik (or m^) as input to fix the physical scale of lattice spacing and the 
u, d and s quark masses. After chiral and continuum extrapolations, the agreement of the calculated 
mass spectrum with experiment is at a 10% level. In comparison with the statistical accuracy of 
1-3% and systematic errors of at most 1.7% we have achieved, this demonstrates a failure of the 
quenched approximation for the hadron spectrum: the hyperfine splitting in the meson sector is too 
small, and in the baryon sector the octet masses and the mass splitting of the decuplet are both 
smaller than experiment. Light quark masses are calculated using two definitions: the conventional 
one and the one based on the axial-vector Ward identity. The two results converge toward the 
continuum limit, yielding m u d = 4.29(14)j^ 79 MeV where the first error is statistical and the 
second one is systematic due to chiral extrapolation. The s quark mass depends on the strange 
hadron mass chosen for input: m s = 113.8(2.3)±| | MeV from m K and m s = 142.3(5.8)lo 2 MeV 
from 77i0, indicating again a failure of the quenched approximation. We obtain the scale of QCD, 
A^-L = 219.5(5.4) MeV with m p used as input. An O(10%) deviation from experiment is observed 
in trie pseudoscalar meson decay constants. 
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I. INTRODUCTION 

Theoretical derivation of the light hadron spectrum 
from the first principles of Quantum Chromodynamics 
(QCD) is a fundamental issue in our understanding of the 
strong interactions. The binding of quarks due to gluons 
cannot be treated perturbatively, and numerical simula- 
tions based on the lattice formulation of QCD, therefore, 
provide a unique means to approach this problem. 

The calculation of the hadron spectrum is made for 
given quark masses and hence it in turn enables us to 
determine the light quark masses, which are the funda- 
mental parameters of QCD. The dynamical scale A of 
QCD is determined by measurements of lattice spacing a 
as a function of the bare coupling constant. Lattice QCD 
also provides us with a method to explore the chiral struc- 
ture which is approximately realized in the real world. A 
further subsidiary verification of QCD may include the 
examination of the decay matrix elements against exper- 
iment. 

Lattice QCD simulations, however, are computation- 
ally demanding, particularly when effects of dynamical 
quarks are to be included. Therefore, since the pioneer- 
ing attempts in 1981 [0J^], the majority of lattice QCD 
simulations have been made within the quenched approx- 
imation in which pair creation and annihilation of sea 
quarks are ignored. In fact such calculations have given 
hadron spectrum in a gross agreement with experiment, 
but clear understanding has not been achieved yet as to 
where this approximation would break down. In order 
to study this point, a calculation with a much higher 
precision is needed. Such a high precision study requires 
accurate controls of a number of systematic errors, which 
is not an easy task even within the quenched approxima- 
tion. The origins of systematic errors include finiteness 
of lattice size, coarseness of lattice spacing and extrapo- 
lations in quark masses from relatively large values. 

The work of GF11 Collaboration carried out in 1991- 
1993 U has advanced the control of systematic errors 
from a finite lattice spacing and a finite lattice size. Tak- 
ing an advantage of large computing power, the GF11 
Collaboration calculated the light hadron spectrum with 
three sets of the coupling constant and three different lat- 
tice sizes at one coupling constant, which is used to take 
continuum limit and estimate finite lattice effects. They 
claimed that the resulting spectrum is in agreement with 
experiment within 6%, the difference for each hadron be- 
ing within their errors. 

We feel that their results need a further verification 
by an independent analysis, since we consider that their 
conclusions depend crucially on the error estimate at sim- 
ulation points, and on a rather long chiral extrapola- 
tion from the region of the pseudoscalar to vector meson 
mass ratio mps/mv = 0.9 - 0.5. Another issue is that 
GF11 simulations were made only for degenerate quarks. 
Masses of strange mesons and decuplet baryons were esti- 
mated using mass formulae, while strange octet baryons 



were not calculated. 

We have embarked on the program to push the calcula- 
tion of the quenched light hadron spectrum beyond that 
of the GF11 Collaboration to answer the posed problems. 
We have aimed to achieve a precision of a few percent for 
statistical errors and to reduce systematic errors to be 
comparable to or smaller than statistical errors. Taking 
the Wilson quark action and the plaquette gluon action, 
simulations are made with lattices of physical spatial size 
L s a w 3 fm for the range of a ~ 0.1 - 0.05 fm. The small- 
est value of mps/mv is lowered to « 0.4. We take an 
advantage of the recent development of quenched chiral 
perturbation theory (Q%PT) which suggests us the 
form of chiral extrapolations. We assume that the light 
u and d quarks are degenerate, but the heavier s quark 
is treated separately giving a different quark mass. 

During the time, the MILC Collaboration carried out 
studies in a similar spirit || using the Kogut-Susskind 
quark action. Because of complications with the spin- 
flavor content of this action, they reported only the nu- 
cleon mass taking and m p as input, leaving aside all 
other hadrons. 

Our calculation was made by the CP-PACS computer, 
a massively parallel computer developed at the Univer- 
sity of Tsukuba completed in September 1996 0. With 
2048 processing nodes, the peak speed of the CP-PACS 
is 614 GFLOPS (614 billion double precision floating- 
point operations per second). Our optimized program 
achieves a sustained speed of 237.5 GFLOPS for the heat- 
bath update of gluon variables and 264.6 GFLOPS for 
the over-relaxation update, and 325.3 GFLOPS for the 
quark propagator solver for which the core part is writ- 
ten in the assembly language Q. The simulations were 
executed from the summer of 1996 to the fall of 1997. 

A brief description of our results has been published 
in Re f. H , and preliminary reports have appeared in 
Refs. (HjJ'Tln this article we present full details of analy- 
ses and results. 

The organization of the paper is as follows: in Sec. || 
the lattice action and simulation parameters are ex- 
plained. In Sec. Ill we present a summary of results 
for the light hadron spectrum, quark masses and meson 
decay constants. Subsequent sections describe details of 
our analyses. In Sec. |y| measurements of hadron masses 
and quark masses at simulation points are discussed. We 
then examine in Sec. [V| the prediction of QxPT for light 
hadron masses against our data. In Sec. VI we describe 
the extrapolation procedure of hadron masses to the chi- 
ral and continuum limits. Compariso ns w ith other stud- 
ies are given in this section. In Sec. VII we discuss de- 



terminations of the light quark masses, and in Sec. VIII 



the QCD A parameter. In Sec. IX results for meson de- 
cay constants are presented. Finally, Sec. [X] presents an 
alternative analysis in which the order of the chiral and 
continuum extrapolations are reversed. Our conclusions 
are given in Sec. 



pa 

Appendices |X] to |Gf 



Technical details are relegated to 
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II. LATTICE ACTION AND SIMULATION 
PARAMETERS 

A. Lattice action 

We generate gauge configurations using the one- 
plaquette gluon action, 



5 9 = f X>Tr(tf P ), 



(1) 



where (3 = 6/g 2 with g the bare gauge coupling constant. 
On gauge configurations, we evaluate quark propagators 
using the Wilson fermion action, 



S q = — il){n)D(n; n, m)ip(m), 



(2) 



+(/+ 7(1 )[/i^ m+ »i, t .}> (3) 

where the hopping parameter k controls the quark mass. 

B. Simulation parameters 

Simulation parameters are summarized in Tables | and 
O. Four values of (3 are chosen so as to cover the range 
of a w 0.1 - 0.05 fm (a' 1 w 2-4 GeV). 

We employ lattices with the physical extent of L s a ~ 3 
fm in the spatial directions. In a previous study, no sig- 
nificant finite lattice size effect was observed for L s a > 2 
fm beyond a statistical error of about 2% p] . For a 
large lattice, the dominant size effect comes from spa- 
tial wrappings of pions whose magnitude decreases as 
Wi s — moo oc exp(— cm^Ls) fll2| ]. For smaller lattices 
squeezing of hadron wave functions enhances the finite 
size effect, leading to a power low behavior, TO£ 3 — Woo oc 
c/Lg [[l3||l4[ |. Assuming the latter behavior, we expect the 
finite-size effects on lattices with L s a rs 3 fm to be about 
0.6%, which is sufficiently small compared with our sta- 
tistical errors. This requires us to use a 64 3 lattice for 
simulations at a w 0.05 fm. 

For the temporal extent of the lattices, we adopt 
L t — (7/4)-Z s . This gives the maximal physical time sep- 
aration of Ltd/2 w 2.5 fm. With our smearing method 
described below, we find that this temporal extent is suf- 
ficient to extract ground state signals in hadron propa- 
gators, suppressing contaminations from excited states. 

For the quark mass, we select five values of k, so that 
they give m PS /m v w 0.75, 0.7, 0.6, 0.5, and 0.4. The two 
heaviest values, which we denote as S\ and S2, are chosen 
to interpolate hadron mass data to the physical point of 
the s quark. The three lighter quarks denoted as Ui, U2 
and U3 are used to extrapolate to the physical point of 
the light u and d quarks, mps/mv = m^/irip = 0.176. 



The quark mass at the smallest value of mps/mv ~ 
0.4 is closer to the chiral limit than that in any previ- 
ous studies with the Wilson quark action, in which cal- 
culations were limited to mps /my>0-5. Reducing the 
quark mass further is not easy. Test runs we carried out 
for nips/my « 0.3 at (3 = 5.9 show that fluctuations 
become too large and the computer time for this point 
alone exceeds the sum of those for the five k down to 
mps/my ~ 0.4. 

Gauge configurations are generated by the 5-hit 
pseudo-heat-bath algorithm |l5[ and an over-relaxation 
algorithm |16|, mixed in the ratio of 1:4. We call the 
combination of one pseudo-heat-bath update sweep fol- 
lowed by four over-relaxation sweeps "iteration". The 
periodic boundary conditions are imposed in all four di- 
rections. The acceptance rate in the pseudo-heat-bath 
step is 82-85% as listed in Table ffl. For vectorization 
and parallelization of the computer program, we adopt 
an even-odd algorithm. 

After 2000-20000 thermalization iterations, we calcu- 
late quark propagators and measure hadronic observable 
on configurations separated by 200 to 2000 iterations de- 
pending on (3, while we measure the gluonic observable, 
such as the plaquette expectation value, at every itera- 
tion. The total number of configurations and their sepa- 
ration are summarized in Table |lj. 

We estimate errors by the jackknife method except oth- 
erwise stated. Tests on the bin size dependence do not 
show the presence of correlations between successive con- 
figurations, and hence we use the unit bin size for error 
analyses. 



Table III shows the number of employed processors of 
the CP-PACS and the execution time required for gen- 
crating and analyzing one configuration. Simulations at 
f3 = 5.9, 6.1 and 6.25 are carried out on subpartitions 
of the CP-PACS computer while at (3 = 6.47 the whole 
system with 2048 processing units is used. 



III. SUMMARY OF RESULTS 

A. Quenched chiral singularity 

Quenched chiral perturbation theory (Q%PT) pre- 
dicts that hadron mass as a function of quark mass m q 
exhibits a characteristic singularity in the chiral limit. 



Data for 



PS 



strongly support the existence of an ex- 



pected singular term 5m q lnm q with 6 ~ 0.1. For vector 
mesons and baryons, the accuracy of mass data and the 
covered range of m q axe not sufficient to establish the 
presence of quenched singularities. 

In Figs. |l| to [|, the QxPT fit is shown by solid lines 
for (a) pseudoscalar meson, (b) vector meson, (c) octet 
baryon and (d) decuplet baryon. The data are consistent 
with the theoretical expectations from QxPT, not o nl; 
for pseudoscalar mesons but also for vector mesons 



>ruy 

0] 
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and baryons |Lg|. We therefore adopt functional forms 
based on QxPT for chiral extrapolations for all cases. 



B. Quenched light hadron spectrum 



We take experimental values of m v = 0.1350 GeV and 
m.p = 0.7684 GeV as input for the mean u, d quark 
mass m Ui d and the lattice spacing a. We use either 
m K = 0.4977 GeV or m^ = 1.0194 GeV for the strange 
quark mass m s . As shown by solid lines in Fig. |5|, hadron 
masses determined at each f3 are well described by a lin- 
ear function of a. 

The quenched hadron spectrum in the continuum limit 
is compared in Fig. |^ with experiment shown by horizon- 



tal bars, with the numerical values given in Table IV. 
Filled symbols use mx as input, and open ones em- 
ploy 77i0. The two error bars show both statistical er- 
ror and the sum of statistical and systematic errors (see 
Sec. VI). Statistical errors are 1-2% for mesons and 2- 
3% for baryons. Estimated systematic errors are at worst 
1.8a of statistical ones, which add only extra 1.7% to sta- 
tistical ones. 

Figure || shows that quenched QCD reproduces the 
global pattern of the light hadron spectrum reasonably 
well, but at the same time systematic deviations exist 
between the quenched spectrum and experiment. An 
important manifestation of this discrepancy is that the 
quenched prediction depends sizable on the choice of par- 
ticle (K or <j>) to hx m s . While an overall agreement in 
the baryon sector is better if is employed as input, 
mx disagrees by 11% (6a), which is the largest difference 
between our result and experiment. 

In the meson sector, the discrepancy is seen in the 
hyperfine splitting, that is too small compared to exper- 
iment. If one uses the mx as input, the vector meson 
masses mx* and m^ are smaller by 4% (4er) and 6% (5a). 
If m<f, is employed instead, mx* agrees with experiment 
within 0.8% (2a) but mx is larger by 11% (6a). 

The smallness of the hyperhne splitting is observed in 
a different way in Fig. [?], which plots m\ — m 2 p S as a func- 
tion of fripg. The hgure shows an approximate scaling 
over the four values of (3. The convergence of data toward 
the experimental point corresponding to (m w ,m p ) is due 
to our choice of these particles as input. Toward heavier 
quark masses, the mass square difference decreases faster 
than experiment, and is about 10% smaller at the point 
corresponding to (mx,mK*) mesons. 

A faster decrease of to? a — m 2 Da can be quantified 



through the J parameter 



PS 

defined by 



J = my 



dmy 
dm PS 



(4) 



A large negative value of the slope seen in Fig. fj] trans- 
lates into a small J as shown in Fig. @; we obtain 



in the continuum limit, to be compared with the experi- 
mental value ~ 0.48 at my /nips — 1.8. 

In the octet baryon sector, the masses are all smaller 
compared to experiment. The nucleon mass is lower than 
experiment by 7% (2.5a). The strange octet baryons are 
lighter by 6-9% with mx as input, and by 2-5% even 
with 777,0 as input. The S-A hyperfine splitting is larger 
by 30% (50%) with mx (m<j>) input, though the deviation 
of 0.8tr (2.3<r) is statistically marginal. The Gell-Mann- 
Okubo (GMO) relation 



-(m N + m s ) = i(3m A + m s ) 



(0) 



J = 0.346(23) 



(5) 



based on first-order flavor SU(3) breaking is well sat- 
isfied, at 1% in both mx and input, though the 
two sides take values (1.04(2) GeV for the m^-input and 
1.09(1) GeV for the m^-input) smaller than experiment 
(1.13 GeV). 

For decuplet baryons, the mass of A turns out to 
be consistent with experiment within statistical error of 
2.0% (0.7a). An equal spacing rule is well satisfied, the 
three spacings mutually agreeing within statistical errors. 
However, the mass splitting is smaller by 30% on average 
compared to experiment for mx input, and by 10% for 
77i0 input. 

The results discussed above are based on QxPT chiral 
fits. In order to see effects of choosing different chiral 
fit functions, we repeat the procedure using low-order 
polynomials in m q , as was done in traditional analyses. 
Chiral fits and continuum extrapolations for this case are 
illustrated in Figs, [j] to || by dashed lines and in Fig. [| by 
open symbols and dashed lines, respectively. QxPT and 
polynomial fits lead to masses which agree within 1.5% 
or 1.6a. The pattern of the quenched spectrum remains 
the same even if one adopts the polynomial chiral fits. 



C. Reversibility of order of the chiral and continuum 
extrapolations 

In order to obtain the physical hadron mass, one con- 
ventionally carries out chiral extrapolation first and then 
takes the continuum extrapolation (we refer to this as 
method A). These two limiting operations can in prin- 
ciple be reversed, and the resulting spectrum should be 
unchanged. An advantage with the reversed limiting pro- 
cedure (method B) is that one need not worry about pos- 
sible 0(a) terms that are present in QxPT formulae at 
finite lattice spacings. 

The light hadron spectrum from the two methods are 
compared in Fig. ^| for the case of the mx input. The 
prediction from the method B denoted by open symbols 
is in good agreement with that of the method A plotted 
with filled symbols within 1.5o\ 

An additional advantage of method B is that the 
hadron mass formula can be obtained as a function of 
an arbitrary quark mass, as shown in the Edinburgh plot 
in Fig. Eel 
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D. Fundamental parameters of QCD 



(0\Vi\V) = e % F v my, 



(14) 



The scale parameter A is the fundamental parameter 
of QCD. We evaluate it in the MS scheme to be 



A^l = 219.5(5.4) MeV, 

MS v ' 



(7) 



when the scale is fixed by m p . 

The definition of quark mass for the Wilson quark ac- 
tion is not unique, because chiral symmetry is broken by 
terms of O(a). We analyze quark masses from two defini- 
tions, the conventional one through the hopping param- 
eter, which we call the VWI quark mass (see Sec. IV C), 
and another defined in terms of the Ward identity for 
axial- vector currents (AWI). 

Figures |ll] and 12 show m uc i and m s renormalized in 
the MS scheme at fj, = 2 GeV as functions of a. The VWI 
and AWI quark masses, differing at finite a, extrapolate 
to a universal value in the continuum limit, in accordance 
with a theoretical expectation. 

A combined linear extrapolation assuming a unique 
value in the continuum limit yields 



m ud 



= 4.29(14)1^ MeV 
m s = 113.8(2.3)1^ MeV 
= 142.3(5. 8)tf-° MeV 



rriK— input, 
m^— input. 



(8) 
(9) 
(10) 



We indicate systematic error arising mainly from chiral 
extrapolations. The value of m s differs by about 20% 
depending on thk or m 4> used as input. The difference 
arises from the small value of meson hyperfine splitting 
in the simulation. 



E. Meson decay constants 



The pseudoscalar meson decay constant fps is defined 



by 



(0|A M |PS) = i Pli fps, 



(11) 



in the continuum notation, where Ap is the axial-vector 
current. The experimental value for tt is / ff — 132 MeV. 
Data for fps are shown in Fig. 13 as a function of a. We 
obtain for physical values 



U = 120.0(5.7) MeV, 

f K = 138.8(4.4) MeV m*- -input. 



(12) 
(13) 



These values are smaller than experiment by 9% (2a) and 
13% (5a), respectively. QxPT predicts that the ratio 
fx I f-K ~ 1 m the quenched QCD is smaller than experi- 
ment by about 30%. This quantity is shown in Fig. 
We obtain fx/ fir — 1 = 0.156(29) which is smaller than 
experiment by 26% (1.9a) as QxPT predicts. 

The vector meson decay constant Fy in the continuum 
theory is defined by 



where and my are the polarization vector and the mass 
of the vector meson V. This is related with another con- 
ventional definition by fy 1 = Fy /my. The experimental 
value of F p is 220(5) MeV, where the charge factor is re- 
moved. Figure |lq summarizes the vector meson decay 
constants. We obtain 



F p = 205.7(6.6) MeV, 

F = 229.4(5.7) MeV m^-input. 



(15) 
(16) 



These values are slightly smaller than experiment; by 
6.7% (2.2a) for F p and by 3.8% (1.6c) for F<j,. 

We summarize meson decay constants in Table [y| 



IV. MEASUREMENTS OF HADRON MASSES 
AND QUARK MASSES 

A. Quark propagators 

We calculate the quark propagator G(m) at a value of 
k by solving 



D(k, n, m)G(m) 



S(n), 



(17) 



where D(k, n, m) is the quark matrix defined in Eq. (^), 
and S(n) is the quark source. In order to enhance the 
ground state signal in the hadronic measurements, we use 
smeared quark sources. For this purpose, we fix gauge 
configurations to the Coulomb gauge as described in Ap- 
pendix 

For the smeared source, we employ an exponential form 
given by 



S(n) 



Aexp( 
1.0 



-B\n\) 



for 
for 



n^0 
n = Q 



(18) 



as motivated by the pion wave function measured by the 
JLQCD Collaboration ]2l[ ]. The smearing radius is ap- 
proximately constant, a/B ps 0.33 fm, over the range 
of (3 we simulate. The quark propagator solver and the 
smearing function are discussed in Appendix M. 



B. Hadron masses 

From quark propagators, we construct hadron propa- 
gators corresponding to degenerate combinations, // and 
/// (/ = s i; s 2, v,i, U2, U3), as well as non-degenerate 
combinations of the type SiUj for mesons, and SiSiUj and 
SiUjUj for baryons; two quarks in baryons are taken to be 
degenerate. We study pseudoscalar and vector mesons, 
and spin 1/2 octet and spin 3/2 decuplet baryons. The 
hadron operators are summarized in Appendix |^. 

Hadron propagators are calculated for all possible com- 
binations of point and smeared sources. At the sink we 
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use only point operators. Effective masses m e ff(t) for 
various combination of quark sources are compared in 
Fig. E6L With our choice of smearing function, m e //(t), 
in almost all cases reaches a plateau from above, suggest- 
ing that the smearing radius is smaller than the actual 
spread of hadron wave functions. The onset of plateau 
is the earliest when the smeared source is used for all 
quarks, and statistical error is the smallest for this case. 
In light of this advantage, we extract masses from hadron 
propagators with all quark sources smeared. 

In order to illustrate the qualityof data, typical effec- 
tive masses are shown in Figs. [17] gOjfor degenerate octet 
baryons at the four (3 values. We extract the ground state 
masses using a single hyperbolic cosine fit for mesons and 
a single exponential fit for baryons, taking account of cor- 
relations among different time slices. In Figs. |l7j-|2C|, the 
horizontal lines are the fit and error, with the range of 
the lines representing the fit range. 

The fit range [t m in,t max ] is chosen based on the fol- 
lowing observations: (1) Value of x 2 /N^ decreases as 
tmin increases and becomes almost constant at a time 
slice which we denote as t x . t x in general depends on 
quark masses. (2) The effective mass shows a plateau 
for t>t x . (3) When t max > t x + 3, x 2 /Ndf is insensitive 
to the choice of t rnax . From these findings, we may use 
t x for t m i n . However, in order to avoid subjectivity in 
the identification of t x and plateau, we adopt i m j„ which 
satisfies the following conditions: 

(1) tmin is larger than t x . 

(2) For each kind of particle at each /3, t m i n is common 
to all quark masses. 

(3) For each particle, values of t m i n in physical units is 
approximately constant for all (3. 

We find that these conditions are satisfied by t m i n sjI.O 
fm, 0.7 fm and 0.5 fm for pseudoscalar mesons, vector 
mesons and baryons, respectively. 

The largest time slices t max for vector mesons and 
baryons at j3 = 5.9, 6.1 and 6.25 are chosen by the 
requirement that the error of propagator at t max does 
not exceed 3%. We employ the same criterion for vector 
mesons at j3 = 6.47. For baryons at j3 = 6.47, the fit- 
ting interval becomes too narrow for light quark masses 
if we employ the cut at 3%. We therefore adopt the cuts 
at 3.2% for octet baryons and 4% for decuplet baryons, 
respectively. 

Values of t max are determined with a different strat- 
egy for pseudoscalar mesons, for which we make chiral 
extrapolations taking account of correlations among dif- 
ferent quark masses. (Correlation among different quark 
masses is ignored for other hadrons.) A large value of 
traax results in full covariance matrices with too large 
dimensions. Such matrices frequently have quite small 
eigenvalues due to statistical fluctuations, and lead to a 
failure of the convergence of the fit. In order to avoid in- 
stability of chiral extrapolations, we determine t max by 



trial and error. We adopt t max = 28, 35, 25 and 35 for 
(3 = 5.9, 6.1, 6.25 and 6.47. 

All hadron masses are stable under a variation of the 
fit range. As an example, the fits with the range [t m j n + 
2, t max ] give results consistent within la. Uncorrelated 
fits yield masses consistent with those from correlated 
fits within la for most cases, although the differences are 
about 2a for some cases. 

Errors are determined by a unit increase of x 2 ■ Jack- 
knife errors are similar in magnitude, the difference being 
at most about 25%. The values of \ 2 l^df turn out to be 
consistent or smaller than unity within errors estimated 
by the jackknifc method, meaning that our one-mass fits 
reproduce the hadron propagator data well. This sug- 
gests that the contamination of excited states are well 
suppressed in our fits. Our results for the hadron ma sses 
are repro duced in Tables [v| for mesons and Tables VII 
and VIII for octet and decuplet baryons. 



C. Quark masses 

The bare quark mass is conventionally defined by 
1/1 1 



rn 



VWI(0) 



2 \ Hi f"\/f J 



(19) 



where k c is the critical hopping parameter at which the 
pion mass vanishes. This quark mass is called the vector 
Ward identity(VWI) quark mass, since the divergence of 
the vector current is proportional to the mass difference 
of quark flavors in the current, which looks similar to 
Eq. (@). 

The definition of the AWI quark mass is based on the 
Ward identity for axial- vector currents |22] ]. For the fla- 
vor combination (/, g) , it takes the form, 



(VMn)0) = (m f + m g )(P(n)0) 
+ (SO)+0(a), 



(20) 



where 
.4 



is the axial- vector current and 

P ( n ) = InlSdn 



(21) 



(22) 



is the pseudoscalar density. In Eq.(po|), V M F(n) = 
F(n)—F(n—fi) is the backward lattice derivative, and SO 
is the response of the operator O under the chiral trans- 
formation. The 0{a) term is due to explicit violation of 
chiral symmetry with the Wilson quark action. 



In order to extract 



AW 1(0) 



we use the relation EJ 



AWI(0) 



m f -<+m^W=lirn 



(V 4 A 4 (t)P(0)) 
(P(t)P(O)) 



(23) 
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where 



Ai(t) = Y,Mn={n,t)), 

n 

P(i) = J>(n = (fM)). 



(24) 
(25) 



are projected to zero spatial momentum. The right-hand 
side of Eq. (E3h is evaluated by a constant fit to the ratio 



(V 4 A 4 (t)P(0)) 



sym 



(P(i)P(O)) 



(26) 



where the suffix sym implies a symmctrization of the 
derivative defined by 

(V 4 A 4 (t)P(0)) 8yin = 

(G(t) - G(t - 1) + G(L t -t)- G(L t - t - l))/4 (27) 

with G(t) = (A4(t)P(0)) . For the pseudoscalar oper- 
ator at the origin, we use the smeared source for two 
quarks. Figures |l] and ^ illustrate our data for Eq. j2^ ) 
at (3 = 5.9 and 6.47. The constant fit is carried out with- 
out taking account of correlations between different time 
slices or between the two correlators. Fit ranges are de- 
termined from the plateau of the effective mass. The 



results for 



blc IX 



AWI(0) 



AWI(0) 



are summarized in Ta- 



V. QUENCHED CHIRAL SINGULARITIES 

The chiral extrapolation is conventionally carried out 
assuming a low-order polynomial in quark masses. Chi- 
ral perturbation theory |24], however, predicts a singular 
quark mass dependence in the chiral limit due to the 
presence of massless pions. The singularity is expected 
to be enhanced in quenched QCD |^,|| since the r\i meson 
is also massless in this approximation. In order to choose 
the functional form for the chiral extrapolation, we ex- 
amine whether hadron mass data are consistent with the 
predictions of QxPT. 



where terms proportional to (mi — 1712) 2 are absent [ p5[ . 
The logarithmic term proportional to 6 represents the 
leading quenched singularity. To the leading order in the 
1/N C expansion in terms of the number of colors N c , 6 
is related to the pseudoscalar meson mass and the pion 
decay constant / by 



V 



K 



24tt 2 / 2 



(29) 



Taking the experimental values of / and pseudoscalar 
meson masses, one finds 8 ss 0.2 as a phenomenological 
estimate. The constant a$ in Eq.(^) represents the co- 
efficient of the kinetic term of the flavor singlet meson 
field, which is sub- leading in terms of 1/N C . The mass 
formula also contains a scale A x of O(l) GeV. The pa- 
rameters such as (5, a$ and A x may differ in quenched 
QCD from those in the full theory. 

In order to see whether pseudoscalar mass data exhibit 
the presence of the logarithmic terms, we investigate the 
ratio m 2 PS 12 /(mi + m2) as a function of 7711+TO2. We use 
the AWI quark mass rather than the VWI mass to avoid 
uncertainties due to the necessity of choosing re c , which in 
turn depends on details of chiral extrapolations. Another 
important point with the use of the AWI quark mass 
is that it is free from chiral singularities which cancel 
between the numerator and the denominator in Eq. ( [23] ) 

In Fig. we plot m 2 PS l2 /{ra\ + 7712) as a function of 
mi + m-2- The AWI quark mass is converted to renor- 
mali zed v alues in the MS scheme at the scale 2 GeV (see 
Sec. VII) and the ratio is translated to physical units. 
This enables us to compare the results at different values 
of (3 with the same scale of the figure. We find a clear 
increase of the ratio towards the chiral limit at all values 
of (3 as expected from Eq. ( p^ ) . 

In order to make a more quantitative analysis, we con- 
sider the ratio defined by 



y 



2mi 



PS, 12 



2m 2 



'PS.12 



777,1 + 7772 m ps 11 m l + m 2 



(30) 



PS,22 



A. Mass ratio test for pseudoscalar meson 

For pseudoscalar mesons made of quarks with masses 
mi and 777,2, QxPT predicts the mass formula j^,^) given 
by 



rn 



PS, 12 



A(mi + 777 2 )|f -£[ln— — + 



ImiA 7772 "72 n 

1 In 

7772 — mi nil 



a$A r 2m\A 2m2A 

7771 m 7^5 h 777,2 In ■ 



12tt 2 / 2 

77717772 . 777 2 - 

In 



A 2 

v 



A 2 

X 



7772 — m i mi ) 

+B (mi +m 2 ) 2 + 0(m 3 ,<5 2 ), 



(28) 



Assuming that 6 and a$ as well as quark masses are 
small, we expect 



y = l + Sx + a x z + 0(m 2 ,(5 2 ), 



where 



a^A 2 
12tt 2 / 2 



and the parameters 



7771 + 777 2 . / "7,2 \ 

x = 2 H In — . 

mi — 7772 \mi / 



and 



(31) 



(32) 



(33) 



f / 277717772 777 2 , 

Z = — In 777l - 777 2 , (34) 

A \7772 — 7771 mi 



7 



0{m? q In m q 



represent terms of 0(m q In m q ) and 
Eq. (§). 

We plot y as a function of a; in Fig. E4l with numerical 
values listed in Table |x|. The points fall within a narrow 
ridge limited by two lines y rs 1 + (0.08-0.12) a;. A one- 
parameter fit ignoring ax z and higher order terms yields 
S = 0.10-0.12 depending on (3 as listed in Table [xf 

A two-parameter fit keeping the axz term requires the 
value of A. We estimate A ss 3 GeV for all (3 values from 
data at mps /my— 0.75 and 0.7, assuming m PS = 2Am q . 
Setting / = 132 MeV, we obtain 5 and a$ given in the 
right column of Table XI. For the two finer lattices at 
/3 = 6.25 and 6.47, a$ is consistent with zero and 8 ~ 0.1. 
At the coarser lattices the values of a$ and 5 are not 
stable. 

Further data are needed to pin down precise values of 5 
and a$. We consider results at finer lattices, closer to the 
continuum limit, are more reliable, and take 5 — 0.10(2) 
and a$ = as our best estimate. 

The ratio test for the existence of quenched logarithm 
terms was originally proposed in Ref. [^), in which one 
plots yo defined by 



yo 



PS, 12 



2mi 



m PS 11 mi + 7712 



(35) 



C. Pseudoscalar meson mass fit 

The parameter 5 is estimated also by fitting pseu- 
doscalar meson mass data to Eq. (|2"s| ) assuming a* = 
for all (3. The AWI quark mass introduces errors in the 
fit variable. Therefore the VWI quark mass, rrhj WI ^ = 
(1/k— 1 / k c )/2, is employed taking n c as a parameter. We 
carry out fully-correlated fits described in Appendix 
independently for degenerate and non-degenerate cases. 

A noticeable property of the QxPT formula Eq. ( psj ) 
is that A, S and A x cannot be determined simultane- 
ously because the three conditions to minimize \ 2 are 
not mathematically independent. A possible method is 
to fix /a = 2 A 2 /A 2 . Values of A and S depend on the 
choice of /a, while n c and B, as well as % 2 and the fit 
curve are independent. We consider f\ = 4, 8, 16, and 
32, which correspond to A x » 1.32, 1.00, 0.76, and 0.57 
GeV, respectively.^] This range of A x contains a natural 
scale for chiral perturbation theory, A x = m , or 1 GeV. 



as a function of xq defined by 



The results are summarized in Table XIII. The value 
of S is stable against a variation of /a and (3, and is 
consistent within 2a with our estimate 0.10(2) from the 
mass ratio test. 



D. Comparison with other results for S 



Xq 



m 2 



mi — !7i2 



In 



(36) 



The relation yo = 1 + 8 • follows if we ignore 0(m q ) 
term in Eq.([2^). As shown in Fig. |25|, however, yo sys- 
tematically varies with quark masses, suggesting a con- 
tribution from the 0(m q ) term. The double ratio defined 
by Eq. (30) is designed to cancel the 0(m q ) terms, and 
hence is more effective to observe the quenched singular- 
ity. 



The value of S has recently been estimated by other 
groups. The FNAL group reported S = 0.065(13) g| 
using the clover quark action, and the QCDSF Collabora- 
tion obtained S w 0.14(2) |29| with a non-perturbatively 
improved clover action. These estimates are consistent 
with ours. 

It has been pointed out in Ref. that the x-y corre- 
lation seen in Fig. [24] may be reproduced with a small 5 « 
0.03 - 0.07, if a$ ~ 0.5. Our data for large (3, however, 
do not seem to be compatible with such a large a$. 



B. Ratio test of pseudoscalar meson decay constant 



E. Vector meson and baryon masses 



Quenched chiral singularities are also expected in me- 
son decay constants p7| . Let ff g be the pseudoscalar 
meson decay constant for the flavor combination (/, g) . 
In Ref. §711, the ratio 



Vf 



fl2 



/ll/22 

was shown to satisfy the relation 



l- 6 -x 
2 



(37) 



(38) 



Vf : 

with x defined by Eq.(3 
of yj are listed in Table 
is summarized in Fig. 
agreement with the value from the mass ratio analysis. 



20 



where a$ is set to zero. Values 
and the ratio test with Eq. j37| ) 
We find S w 0.08-0.16, in 



QxPT predicts singularities of the form 0(mps) ~ 
0(^/m q ) for vector mesons and baryons [ fl7|JT8]l . The ra- 
tio tests similar to those for the pseudoscalar mesons in- 
dicate that the coefficient of the 0(mps) term is non- 
vanishing both for vector mesons and baryons. It is dif- 
ficult, however, to reliably estimate the coefficients from 
the ratios because of large errors. Direct fits of mass data 
to QxPT formula are also difficult as they are not very 



1 These values of A x in physical units are computed using A 
determined by a degenerate fit at j3 = 5.90. We confirm that 
13 dependence of A is very weak if translated into physical 
units and that degenerate and non-d egen erate fits lead to A 
consistent with each other. See. Table PCIB. 



8 



stable. While our data are consistent with QxPT, statis- 
tics and the range of the quark mass in our study do not 
allow conclusive results. Our tests of the QxPT mass 
formulae for these cases are described in Appendix |]. 



VI. HADRON MASS SPECTRUM 
A. Chiral fits 

Chiral fits of t he p seudoscalar meson mass are already 
described in Sec. VC , and are shown in Fig s. 0-(a) - ||-(a) 
with parameters summarized in Table KnL Comparisons 
of various fit functions for pseudoscalar meson masses are 
given in Appendix ||. 

For vector mesons and baryons, we choose the pseu- 
doscalar meson mass as the variable to represent the 
quark mass dependence. For vector mesons we adopt 



m v ,i2 = m v 



r 3. . - m?i , 

+^rU m " + m22 + 2 ^ ri 

6 2 m.22 — 



C 

+ ^-( m ll + m 22)> 



(39) 



where m/ s is the pseudoscalar meson mass with the 
quark flavor combination (f,g). The coefficient C±/2 is 
proportional to S while the C\ term is present in ordinary 
%PT. For octet baryons, we employ JTa] 



m s = rrio 



1 



{AF 2 w uu - 4(D - F)Fw u d + (D - F) 2 w ss } 



-Ab F m 2 uu + 2(b D - b F )m 2 ss 
+(2£> 2 /3 - 2F 2 )v au 
+{2D 2 /3-ADF + 2F 2 )v ud , 



(40) 



1 



|{(4£>/3 - 2F) 2 w uu + {D/3 + F) 2 w ss 

-2(4,0/3 - 2F)(D/Z + F)w ud ] 
+4(26^/3 - b F )m 2 uu - 2(b D /S + b F )m 2 s 
+ {2D 2 /9 - 8DF/3 + 2F 2 )v uu 
+{10D 2 /9-4DF/3-2F 2 )v ud , (41) 

for S-type and A-type cases. For decuplet, the formula 
reads, 

5H 2 

m D = m° D + — —{Awuu + 4w„d + w ss ) 
loz 

C 2 



+c(2m 2 uu + m 2 J. 



(42) 



w fi 



v fg 



-27n5- 



l ff 



" L fg 

87T/2 



(43) 
(44) 



Here 



The 0{m PS ) terms are not included in Eq. ( |39| ) for 
vector mesons and in Eq. ( p2[ ) for decuplet baryons. The 
octet-decuplet coupling terms are also ignored in Eqs.(|4fj|) 
and (^lj) for octet baryons. These choices are made be- 
cause fitting parameters are not well determined if these 
terms are introduced (see Appendix |f]), and the dropped 
terms have small effects for the spectrum. Fittings with 
and without them are compared in Fig.^7] for degenerate 
masses at j3 = 5.90. The two types of fittings reproduce 
data equally well. The difference remains small at the 
physical point, at most 5% (5a) at finite lattice spacings 
and at most 1.2% (1.3cr) after the continuum extrapola- 
tion. 

We set 8 = 0.1 and a$ = as suggested from the 
pseudoscalar case. These choices do not affect the fits 
for vector mesons and decuplet baryons: a non- vanishing 
a<s> leads to 0(m^ g ) effect which is not included in the fit 
function, and a change of 6 is absorbed by a redefinitions 
of the parameters. 

For the nucleoli mass, dropping 0{m PS ) terms in 
Eq. ( pi(]| ) would lead to a positive curvature (concave func- 
tion), which contradicts the data that show a negative 
curvature (convex function) in Figs. 0-(c) - §-(c). There- 
fore we include 0(m PS ) terms for octet baryons. The 
coefficient of 0(m PS ) terms is affected by the choice of 
S and a$. We study the effect of the uncertainty of 5 
and q$ on the resulting octet masses by varying S from 
0.08 to 0.12 and a$ from -0.7 to +0.7. The change of 
(5 in this range results in 0.4% (1.3er) difference at finite 
lattice spacings and 0.3% (0.3ct) in the continuum limit. 
The change of a$ leads to differences of 2.9%(4.7<r) and 
2.2%(1.4ct), respectively. We also fix / = f„ = 132 MeV. 
Changing / to fx = 226 MeV affects octet baryon masses 
by at most 2.5a at finite lattice spacings, and 0.5a in the 
continuum limit. Artifacts of fixing these parameters are 
sufficiently small at least in the continuum limit. 

Fits are made to degenerate and non-degenerate data 
together. Because the size of covariance matrix becomes 
too large and the matrix elements cannot be determined 
reliably, we do not include correlations among different 
quark masses. 

Fits for vector mesons and baryons are shown in Figs. |l| 
to ||. Parameters are summarized in Table pCIV for v ector 
mesons, m Table [xy| for octet baryons and in Table |XVI| 
for decuplet baryons. As we see in the figures and tables, 
QxPT fits reproduce data at all (3. 



B. Hadron masses at the physical point 

The extrapolation and interpolation to the physical 
point is made as follows. For the case of the m,K- 
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input, we determine mps[n u d, si) and mps(K u di S2) from 
non-degenerate fits to pseudoscalar mesons, and interpo- 
late them linearly in terms of the s quark mass so that 
mps(, K ud, K s) takes the experimental value of tuk- For 
the m^-input, we first determine the mass m ss of the 
degenerate pseudoscalar meson consisting of two strange 
quarks from the vector meson mass fit and evaluate k s . 
We then make a linear interpolation of rnps(n u d, si) and 
mps{n u d, s 2 ) to find m K = mps(K u d, K s ) = m us . Values 
of m uu , m ss and m us then lead to the predictions for 

other hadron masses. 

Hadron masses in lattice units are listed in Table XVII . 
We include results for fictitious hadrons such as u rj s = ss" 
(pseudoscalar meson consisting of two strange quarks), 
A ss (A-like baryon consisting of two quarks with m s and a 
light quark with m u d), and N sss (octet baryon consisting 
of three quarks with m s ). Hadron masses translated to 



physical units are compiled in Table XVIII 



C. Continuum extrapolation 



For our lattice actions scaling violation is given by 



m(a) = m(0)[l + Sa + (S'af + 0(a 3 )]. 



(45) 



Hadron mass data in Fig. || are fitted well without the 
0(a 2 ) and higher order terms (x 2 /Ndf < 1.6). Hadron 
masses in the continuum limit are give n in Table IV. 
Statistical errors are about 1-3%. Table |XT| lists S of 
the term linear in a. The S for baryons are larger than 
those for mesons. Lighter baryons have larger values of S 
for both octet and decuplet. The nucleon has the largest 
S of about 280 MeV, with which the scaling violation Sa 
is about 10% in the middle of our range of lattice spacing, 
a = 0.075 fm. 

A fit retaining a quadratic term leads to S and S' ill 
determined with the magnitude of errors comparable to 
the central values. The masses in the continuum limit 
have large errors of 13% which are 5 times those from 
the linear extrapolation. Within these errors, the con- 
tinuum results from the quadratic fit are consistent with 
those of the linear one. With only four points of lattice 
spacings, we are not able to test effects of higher order 
terms further. 

We therefore estimate the systematic errors from the 
0{a 2 ) terms by an order estimate, assuming S' ~ S and 
substituting a — 0.075 fm as abovej^] The errors esti- 
mated in this way nor mal ized by the central values are 
summarized in Table XX under the column u O(a 2 ) er- 
ror" . We find their magnitude to be quite small. Even for 
the nucleon with the largest scaling violation, the 0(a 2 ) 



2 For an estimate of 5", we fit deviations of m(o) from the 
linear fit using a pure quadratic function of a. This method 
gives 5" of O(50 MeV) with error of O(100 MeV). 



error is about 1%. Thus, unless S' is unduly large, 0(a 2 ) 
systematic errors would not exceed a percent level. This 
is much smaller than the deviation between the calcu- 
lated quenched spectrum and experiment. 



D. Results from polynomial chiral fits 

Polynomial chiral fits are carried out to degenerate and 
non-degenerate data separately, fully incorporating the 
correlation among different quark masses. We employ 
quadratic polynomials in terms of the VWI quark mass, 
except a cubic polynomial for degenerate octet baryons. 
The fitting procedure is described in Appendix The 
fits are plotted in Figs. |] - f| by dashed lines. Hadr on 
masses at the physical points are listed in Table XXI. 

Extrapolating the results to the continuum limit, 
hadron masses from polynomial chiral fits are also fit- 
ted well by linear functions in a with x 2 /Ndf< 1.7. The 
continuum extrapolation is shown in Fig. M by dashed 
lines. Masses in the continuum limit are given in Table 



XXH 



At the four /3 values, the difference in hadron masses at 
the physical point between the QxPT fit and the poly- 
nomial fit is at most 3%. In the continuum limit the 
differences are within 1.5% (1.6 cr), as listed in the col- 
umn "chiral fit error" in Table XX. This difference is 



sufficiently small so that it does not alter the pattern 
of deviation between the quenched spectrum calculated 
with Q^PT chiral extrapolation and the experimental 
spectrum as shown in Fig. E3. 



E. Systematic error and final results 

Total systematic error for the mass spectrum is esti- 
mated by adding in quadrature the error from continuum 
extrapolation ( u O(a 2 ^ 



error" in Table XX), that from 
chiral extrapolations ("chiral fit error" in the table), and 
a 0.6% error for finite size effects. 

Our final results for the quenched light hadron spec- 
trum including the systematic error are summarized in 
Fig. I and Table [fV. 



F. Comparison with previous results 

1. Meson hyperfine splitting 

The GF11 Collaboration || calculated hadron masses 
with the m^-input using lattices with L s a ~ 2.3 fm. 
The chiral and continuum extrapolations are made with 
a linear form. Based on a finite size study at p = 5.7 
with L s a « 2.3 fm (L s = 16) and w 3.4 fm (L s = 24), 
they corrected the continuum results for finite size effects. 
They claimed that the hyperfine splitting between K and 
K* is consistent with experiment. 
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This differs from our small hyperfme splitting. We 
compare our data (filled symbols) and the GF11 data 
(open symbols), both with to if -input, in Fig. For 
the GF11 data, the results from the larger lattice with 
L s a rs 3.4 fm are also shown (open squares) at a ~ 0.7 
GeV -1 , and the continuum estimates before and after 
the the finite size correction are shown at a = 0. 

We observe that all data for K* and <fi at finite a are 
nearly consistent with each other. The difference in the 
continuum limit is due to a steeper slope of the GF11 data 
for the continuum extrapolation, arising from small val- 
ues of rriK* and at (3 = 5.7 on the lattice of L s a ~ 2.3 
fm (the rightmost triangles). If we adopted the data from 
the L s a ~ 3.4 lattice (open squares), we would obtain a 
continuum value in agreement with our result. 

In Ref. ||, the discrepancy between L s a ~ 3.4 and 
2.3 fm is considered as finite size effects. However, since 
data at smaller a are consistent between L s a w 3.0 fm 
(our data) and 2.3 fm (Ref. ||), it is not clear whether we 
can attribute the difference simply to finite size effects. 
The conclusion of the GF11 critically depends on their 
data at [3 = 5.7, for which we suspect an underestimation 
of errors. 



2. Nucleon mass 

In previous calculations at (3 ~ 5.7-6.2 with 
mps/my>0.5, nucleon masses are significantly higher 
than experiment at finite a pjsifl. The GF11 claimed 
agreement with experiment after the continuum extrapo- 
lation and the finite size correction. In the present study, 
however, we find the nucleon mass to be smaller than the 
previous estimates even at finite a. Extrapolating to the 
continuum limit, we find the nucleon mass to be smaller 
than experiment by 7% (2.5<r). See Fig. 29 where our 
data and those of the GF11 are compared. 



The origin of our small nucleon mass at a finite a is the 
negative curvature in 1/k toward small quark masses, as 
observed in Figs. [I] - [I| This trend becomes manifest 
only when the quark mass is reduced to mps/mv ~ 0.4 
while sustaining statistical precisions. In fact, a linear fit 
of our data at mps/mv>0.5 gives a larger nucleon mass 
consistent with the previous results. 



3. Masses of H * and Q 

The GF11 reported the masses of S* and from the 
m^-input higher than experiment by 3-5%. In contrast, 
Fig. ^| shows that our masses are smaller than experiment 
by a similar magnitude. 

The origin of these differences can be seen in the top 
panel in Fig. ||^. While the results from the two groups 
are consistent at a w 0.5 GeV -1 , the continuum extrap- 
olation is different, especially for f2. As in the case of the 
meson hyperfme splitting, the continuum extrapolation 



of the GF11 is critically affected by the data at (3 = 5.7 
and Lm ~ 2.3 fm. 



4- Comparison with staggered quark results 

The negative curvature of the nucleon mass has also 
been reported in Ref. ||, in which the nucleon mass 
for the staggered quark action is calculated down to 
mps/ m v ~ 0.3-0.4. However, our result m^v = 878(25) 
MeV in the continuum limit obtained from the Wilson 
quark action is smaller than = 964(35) MeV ||] from 
the staggered quark action by about 2.5 a. 

It has been pointed out in Ref. |32) that the differ- 
ence in the Wilson and Kogut-Susskind results for the 
nucleon mass exists not only at the physical point but 
even at heavier quark masses for which the discrepancy 
is statistically more significant. In Ref. p2| , the nucleon 
to p mass ratio off the physical quark mass is calculated 
in the continuum limit, using the same method as that 
explained in Sec.[x] below. The ratios for the staggered 
action are larger than ours by about 8% for the whole 
range of the quark mass (see Fig. 4 in Ref. |32]| ). 

The origin of the difference is not explained by finite 
size effects since both calculations employ sufficiently 
large lattices, nor by the chiral extrapolation since the 
difference exists even for heavy quarks. The continuum 
extrapolation is also improbable as the cause, because 
masses calculated at finite (3 are well reproduced by low- 
est order scaling violation of 0(a) for our Wilson results, 
and 0(a 2 ) for the staggered quarks. For the moment the 
origin of the difference is an open issue. 



VII. LIGHT QUARK MASSES 

A. Renormalization factors 

We calculate renormalized quark masses in the MS 
scheme at the scale p = 2 GeV. They are given by 



m y^ = Z m m v q WI ^a-\ 

m AWI = | m */(0) fl -l 

Z P 



(46) 
(47) 



for VWI and AWI masses. The renormalization factors 
Z m [F53), Za 1 34 1, and Z P |35| are estimated with tadpole- 



improved one-loop perturbation theory |36| , by matching 
the lattice scheme to the MS scheme at /1 = 1/a. They 
read 

Z m = 8/t c [l + 0.01a P (l/a)], (48) 
Z A = (l + (0.448 - (47r/12))ap(l/a))/uo, (49) 
Z P = 1- 1.0335a P (l/a), (50) 

with Mo = (Up) 1 / 4 , where (Up) is the plaquette average. 
For ap (1/a), we first compute ap(q*) according to 
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-ln«^» = y«p(?*)[1 



1.18969 a P (q*)}, (51) 



where q* = 3.4018/a, and use the renormalization group 
equation to two-loop order. The running of the quark 
mass from fi = 1/a to 2 GeV is made employing the 
three-loop renormalization group equation [B7J . 



B. Chiral and continuum extrapolations 

For AWI quark masses, we need to carry out chiral 
extrapolation and/or interpolation to the physical point. 
Polynomials in 1/k are used for this since quenched chi- 
ral singularities are absent |25|,^6|]. In fact, as shown in 
Fig. pfl, the ratio of renormalizcd quark masses 



,vwi 



,AWI 



is flat as a function of x = 



,AWI 



,AWI 



,AWI 



(52) 



suggesting 



linear behavior of the AWI quark mass in 1/k. 



A comparison of k c and where 



,AWI 



vanishes 



suggests the presence of QxPT singularity for the pseu- 



-.AWI. 



see 



doscalar meson mass mp$ and its absence for m q 
inset plots of Figs. 0(a) - 0(a). The value of n c from the 
QxPT fit to mps and 1 from a linear fit to mf WI 
agree well, whereas a quadratic fit of mps clearly fails to 
do so. Values for k^ wi and K r . from various fits are com- 



piled in Table XXIII. Numerically, k c from the QxPT 



fit of mp s agrees with n^ WI obtained from a linear or 
quadratic fit with at most 2.8o\ On the other hand, if 
we adopt fits with the quadratic (cubic) extrapolation of 
mp s , the difference between k c and k^ wi increases to 
as much as 17c (12a). 

For actual chiral extrapolation of m^ WI , we employ a 
quadratic fit and enforce n^ WI to agree with k c obtained 
from the QxPT fit of mps, having confirmed their agree- 
ment as described above. This constraint is imposed be- 
cause even a small difference between k^ wi and n c affects 



estimates of : 



,AWI 



at the physical point. We employ 



2mf w/ (°) = B( 9 (1/k - 1/k c ) + B* 9 {1/k - l/n c f. (53) 



Fit ting p arameters Bf 9 and B°i 9 are given in Ta- 



ble XXIV. Fits with various functional forms are hardly 



distinguishable, as shown in Fig. |31j. 

For calculating the s quark mass, we first make a linear 
fit to non-degenerate combination of quark masses, 



l AW m+m A Wm =A £ +A s i/K 



(54) 



Values for m. 



AWI(0) 



and m. 



VWI(0) 



mass point are presented in Table XXV 
transl ated to MS scheme at fi 



at the physical quark 



Quark masses 
2 GeV are given in 
Table [XXVI| , and are plotted in Figs. |ll| and (d| They 
are well reproduced by linear functions in a. The AWI 
and VWI quark masses, which are different at finite a, 
extrapolate to a universal value in the continuum limit, 
as they should. We determine the quark masses in the 
continuum limit by a combined linear extrapolation of 



m 



VWI 



and m£ WI (Table |XXVII| ) 



C. Systematic errors and final results 

To estimate systematic errors from chiral extrapola- 
tions, we consider a quadratic fit to m^ WI taking k^ wi 
as a fit parameter. We then carry out a QxPT fit to 
mps with k c set to the k^ wi , and evaluate VWI and 
AWI quark masses. Chiral fits up to here employ n as an 
independent variable. To evaluate errors in quark masses 
from fits in terms of k, we consider an independent QxPT 
fit to rnpa_as a function of m^ WI without referring to k. 

Figure |32| shows that m u d is sensitively affected by the 
treatment of the chiral limit. At finite a, both m^ WI 
and m^ WI from the alternative fit above shown by tri- 
angles differ from the original ones (circles) far beyond 
statistical errors. Linearly extrapolated to the contin- 
uum, the alternative methods lead to m u d w 4.1 — 4.8 
MeV, depending on the choice of fits. The fit to mps 
as a function of m^ 1 (shown by diamonds) gives the 
lowest value of ~ 3.5 MeV. Taking the maximum of the 
differences between five results and the value obtained 
in the preceding subsection, we estimate the systematic 
error to be +0.51 MeV and -0.79 MeV. 

Systematic errors from chiral fits are not large for m s . 
As Figs. [33] and [34] show, results from various chiral fits 
at finite a agree with each other within at most 3o\ We 
estimate the systematic error in the continuum limit by 
the same method as for m u d- We obtain +5.8 MeV and 
-2.9 MeV for the m K input and +22.0 MeV and -0 
MeV for the input. 

We also investigate uncertainties from various defini- 
tions of the axial- vector current and higher order effects 
in the renormalization factors. For the former, we test 
for the local axial current defined by 

A l ™\n) = f n i 75 ^g n (55) 

with the tadpole improved renormalization factor 

z iocai _ 1 _ .3i6ap(l/a). (56) 



keeping K Si fixed. We then set k = n u d, and calculate 



AWI(0) 



m Awi(o) ^ a ]j near interpolation in terms 
of l/ft Sj . We do not employ a quadratic extrapolation 
since the effect of the quadratic term is negligibly small 
in s-quark mass but increases errors of fitting parameters 
significantly. 



The procedure to calculate the AWI quark mass 
m ^W 1, locality j g dgggj-jbgd m Appendix For the latter, 
we repeat analyses using the MS coupling 



(U P ) 



4s(V«) 



0.134868, 



(57) 
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XXVIII . The very small statistical error of plaquette is 



instead of a pill a). Substituting a^jg(7r/a) into ( p8| ) with [i — ixja yields 

Figures |3q, |3(] and |37] show the data and continuum ex- A^-. Ou r measurements of plaquette are listed in Table 
trapolations. Values of m^ WI (triangles) and m AWi,local 
(squares) are in good agreement. The difference is about 
5% on the coarsest lattice and is smaller on finer lattices 
for all cases of m u d, rn s (m^- input) and m s (m^ input). 
The two values agree in the continuum limit within 1.5cr 
of the statistical error. The results with ap and a-^ 
are compared using filled and open symbols. The small 



difference in 



,vwi 



reflects a small value of one-loop co- 



efficients. On the other hand, the difference in m^ WI 
is about 5% on the coarsest lattice and about 3% on the 
finest lattice, which leads to a difference of 2% in the con- 
tinuum limit, to be compared with the statistical error 
of 2-4%. 

As shown in Figs. |3f| ^ and [57], the central values 
in the continuum limit are contained within the error 
band given by the sum of statistical error and systematic 
one from chiral extrapolations. Therefore we do not add 
the errors from the definition of current and higher order 
effects in the renormalization factors to the estimate of 
the systematic error above. 

Final results are given in Eqs. (^),((|) and (jlC|). We 
note that these numbers are different from those given in 
our earlier publication Q , in which we employed a linear 
chiral extrapolation of m^ WI and corrected for a small 
difference between k c and k^ wi . Values here obtained 
by constrained quadratic chiral fits are our final results. 



VIII. QCD SCALE PARAMETER 

A. Methods and results 

We calculate the QCD scale parameter A^ in the MS 
scheme. In this scheme, the renormalization group co- 
efficients are known to four-loop order. Since the rela- 
tion between the lattice coupling and the MS coupling is 
known only up to two loops [p8| , we employ the expres- 
sion to three-loop order given by 



a ms/a* = (Poy) 



. 01 

20 o exp(- 



1 



)(1 



dl - (3 p 2 



2$ 



-v 



(58) 



«T7s(M)/47r and (3f s 



2857/2 for the MS 



where y - ^ MgV 
scheme of quenched QCD ]39| |. 

We estimate the MS coupling in three ways, using 
the mean- field coupling ckm f pCfl , the plaquette coupling 
ap |36| and the potential coupling ay. 

The mean field (or tadpole improved) coupling cxmf is 
defined by cxmf — &o/(Up), where ag is the bare lattice 
coupling. Using the relation between the MS coupling 
and the bare coupling up to two loop and the per- 
turbative expansion of the plaquette, one obtains 



1 



1 



a MS'( 7r / a ) a MF 



0.30928- 1.95683a A /F- (59) 



ignored in our analyses. 

The potential coupling ay is defined through the static 
qq potential [^6|. Using the plaquette coupling ap de- 
fined by 

- ln([/p) = — ap ( q *)(l - 1.18969a F (g*)), (60) 

with the optimal value of q* = 3.4018/a the po- 

tential coupling is evaluated by 



a v (q*) = a P {q*) + 2.81404af>(g*). 



(61) 



We evolve av{q*) from fi = q* to fj, = it /a using the 
three-loop renormalization group with (3% = 6178.36 
HEl, and calculate 



a Ms(' !T / a ) = Oi V {ir/a) 
-0.82230a y (7r/a) 2 - 2.66504a y (7r/a) 3 



(62) 



We also evaluate a^jg directly from the ap coupling 
using the BLM G3 improved relation 



a m (ir/a) = a P ((J,) + -a P {p) 2 + 0.95465a|>(^) 3 , (63) 
where fi = e 5 / 6 q*. 



Values of A»jg are given in Table XXIX and Fig 



They are fitted well with a function linear in a with small 
X 2 /Ndf = 0.3 - 0.5. The difference among continuum val- 
ues of Aj^jg is smaller than statistical errors. This suggest 
that possible higher order terms in a are not important 
within our accuracy. Since A^g from ay coupling ex- 
hibits the smallest scaling violation, we quote Eq.(0) as 
our best estimate of Aj^g- in quenched QCD. We note that 
the scale is determined from m p . 



B. Comparison with other results 

There are a number of ways to determine the QCD 
scale parameter. The SCRI |43] group obtained A^g- = 
247(16) MeV from the measurement of \/a/A using the 
string tension y/a = 465 MeV estimated from the char- 
monium level splitting. With a recursive finite-size tech- 
nique using the Schrodinger functional the ALPHA Col- 
laboration jl6) obtained A^g = 238(19) MeV, where the 
physical scale is fixed by the Sommer scale |47J r = 0.5 
fm. In what follows we consider A values that would come 
out if we determine the string tension or the Sommer 
scale with the p meson mass obtained in our simulation. 

We borrow a parameterization fill (proposed in 
Ref. (4^]) of \fa and 1/tq obtained with high statistics 
data, and evaluate them at (3 = 5.9, 6.1, 6.25 and 6.47, 
the points of our simulation. Using our m p we evaluate 
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y/a/rrip and l/r^m p as depicted in Fig. 39, where the er- 
rors arising from ^fa and 1/ro are ignored. These values 
show a linear dependence in a, leading to 



(64) 
(65) 



, p = 0.494(16), 
l/(r m p ) = 0.409(13). 



With nip — 768.4 MeV, we obtain in the continuum limit 
= 380(12) MeV and r = 0.628(20) fm, which exhibit 
10—20% deviations from the usually accepted values. 

With these scales the SCRI result, A^g = 247(16) 
MeV, is converted to A^g = 202(13) (7) MeV, while 
A^ = 238(19) MeV, to A^ = 



the ALPHA result, 



, MS - 238(19) MeV 
189(15)(6) MeV. Here, the first errors come from those 
quoted in the original literature, and the second are from 
the error of our p mass measurement. Figure [}8] com- 
pares our estimate of A with those obtained by SCRI 
and ALPHA, and with those we have re-evaluated con- 
sistently using the p mass (labelled as "translated"). It 
seems that there are discrepancies somewhat beyond er- 
rors. Three values of A obtained with the same scale 
(m p in this case) are supposed to be those in the contin- 
uum limit and hence discrepancy, if it exists, should be 
resolved within quenched QCD. 



IX. MESON DECAY CONSTANTS 



A. Pseudo-scalar meson decay constants 



We extract pseudoscalar meson decay constants fps 
from the correlation function (A l £ cal (t)P(0)). The value 



of j ps is related torn, 



f PS a = Z A fpg, 



defined in Eq. (p2j): 

(66) 



f(0) _ 
J PS ~ 



2mpsaCp S 



X (TOy 



{mp S a) 2 

AWI,local(0) 



AWIlo 



22(0) 



(67) 



for the flavor combination (/, g) , where Cp S is the ampli- 
tude of propagator with the point quark source defined 
by 



(P p (t)P p (0))^C P s exp(~m PS t). 



(68) 



To avoid the direct use of point source propagators, 
which are noisier than the corresponding smeared source 
propagators, we apply the following procedure: the am- 
plitude Cp S for the smeared source propagator is already 
obtained from the meson mass fit 



(P p {t)P s (0)) w C| s exp(-mp S at). 



(69) 



The ratio ?y = Cp S /Cp S can be obtained from an addi- 
tional fit 



(p p (t)p p m 

(P p (t)ps(0)) 



(70) 



Note that Eq. (]7^) is the only new fit which is necessary 
to calculate fps- We illustrate a typical fit in Fig. [io[ 

For the renormalization constant Za, we employ the 
tadpole-improved ^,^] one-loop formula [B5| , which is 
factorized into two parts, 



Z A 



rylocal ry 
^ A K 



where 



Z K 



3k/ 
4k c 



% K g 

4k c 



(71) 



(72) 



Table XXX summarizes the results for Z K f p g at the sim- 
ulation points. 

For / = Z K fpg, non-degenerate data are well repro- 
duced by a linear function of l/n Uil where K Ui is the 
hopping parameter for the u, d quark, while degenerate 
data show a slight negative curvature (Fig. ^lj). There- 
fore, we employ a quadratic polynomial function of 1/k 



f = A d 9 + A d 1 9 /K + A^/K 2 



(73) 



for the degenerate case, and a linear function of 1/ K Ui 

f = A s >+A s 1 i / Kui (74) 



for the non-degenerate case (it,, Sj). Table XXXI presents 

Z K f P s at the physical point. 

Multiplyin g Z l £ cal and a -1 , we obtain fp$ in physical 



units (Table XXXII ) . Data are well reproduced by a 
linear function in a, as shown in Fig. [l^. Accordingly, 
we obtain small x 2 /N d f = 0.45, 0.86 and 1.38 for /„., f K 
(rnif-input) and fx (m^-input), respectively. Our final 
results for fps are summarized in Eqs. (^), ( p^|) and 
Table 0. For the ratio f K /U - 1, we obtain 0.156(29) in 
the continuum limit from a linear fit of Fig. |l4| . 



B. Vector meson decay constants 

Vector meson decay constants Fy are extracted from 
the correlation function (Vi(n)Vi(0)) of the local vector 
current 

Vi(Tl) = Jnli9n- (75) 

Since Vi is the vector meson operator, we obtain 

^ (0) (76) 
(77) 



Fya = Zv\/2Kf\/2KgFy , 
4 0) = ^2C p /m v a, 



where C p is the amplitude of vector meson propagator 
with the point source. Employing a method similar to 
that used in the previous subsection, we first obtain the 
amplitude of smeared propagator Cy and then fit the 
ratio rj of the point and smeared propagators to calculate 
C p = r,Cl 
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We adopt a non-perturbative definition |50| given by 

(v c W(o)) 



Zv = lim 



(78) 



f^oo (V(t)V(Q)) ' 

fi n ) = oj/^+A^n./iCTi + l)g n + fnUnAli ~ 1 )5«+a} 



where 



is the conserved vector current. 
Eq. are shown in Fig. E|. 



(79) 

Examples of the fit for 
Values of dep end 
We 



little on the quark mass, as shown in Table XXXIII 
provide Fya at the simulation points in Table XXXIV 



Chiral extrapolations are carried out in a similar man- 
ner to fps- Dat a arc fit by a linear function of 1/k (see 



Table XXXV presents F p and (m^ input) 



in units of GeV. We extrapolate them to the continuum 
limit linearly in a, as shown in Fig. [l5|. Final results for 
Fy using the non-perturbative renormalization constant 
are summarized in Eqs. (p"5|), (16) and Table [v|. 

For comparison, we study the renormalization fac- 
tor estimated by tadpole-improved one-loop perturbation 
theory; 



F\ 



y P a = yjl — 3k//4k c 



1 — 3k 9 /4k c 



x(l - 0.82ap(l/o))F, 



(0) 



(80) 



Continuum extrapolations, plotted in Fig. |15|, show that 
Fy p exhibits scaling violation much larger than Fy . Al- 
though the difference of Fy p and Fy at finite a becomes 
smaller towards the continuum limit, Fy P extrapolates 
to a value slightly smaller than Fy. The coefficient of 
ap in Eq. (pOI) is rather large, when compared to that 
for Z l ^ cal in Eq. ( |56| ) used for fp$. Higher order terms 
may be important to extract Fy p in the continuum limit 
from the region of our lattice spacing. Therefore we take 
Fy determined with the non-perturbative renormaliza- 
tion factor as our best estimate. 



X. LIGHT HADRON SPECTRUM AS 
FUNCTIONS OF QUARK MASSES 

We carry out an additional analysis in which the con- 
tinuum extrapolations are made before the chiral fits, and 
obtain the hadron spectrum in the continuum limit as 
functions of quark masses. A motivation of this analysis 
concerns the question of whether QxPT mass formulae at 
finite a may suffer from lattice artifacts. Because QxPT 
parameters obtained in Sec. VI are smooth in a, we ex- 
pect that 0(a) terms vanish smoothly toward the contin- 
uum limit. The alternative procedure provides us with a 
more direct test of the QxPT formulae in the continuum 
limit. Furthermore, quark mass or mps/my dependence 
of hadron masses in the continuum limit can be used to 
estimate the size of the scaling violation in future calcu- 
lations with improved actions at fixed a, without difficult 
chiral extrapolations. 



A. Continuum extrapolation 

We take the continuum extrapolation of hadron mass 
data at finite quark mass. To do this, we first interpolate 
or slightly extrapolate data to mps/my — 0.75, 0.7, 0.6, 
0.5, and 0.4 at each (3. These values of mps/my are 
close to our raw data points such that the errors and 
uncertainties from fits are small. In practice, we use the 
QxPT formulae adopted in Sec. VI. We also repeat the 



whole procedure using polynomial chiral fits to estimate 
systematics of choosing different fit functions. Errors for 
interpolated/extrapolated values are a factor 1-3 larger 
than those for the raw data. 

We then extrapolate the hadron masses to the con- 
tinuum limit at each mps/my. In order to calculate a 
relative value of a as a function of f3, we use the vector 
meson mass my at mps/my = 0.75, which we denote 

(0.75) 



as m\ 
by to 



(0.75) 



i.e., masses at each mps/my are normalized 



and extrapolated linearly in m 



(0.75), 



to the 



continuum limit. A value of 



(0.75) 



in physical units is 



not necessary at this step. We provide in Tables XXXVI 



XL 



normalized hadron masses at each {3 obtained by 
the QxPT formulae and those in the continuum limit. 
Values of x 2 /Ndf for the continuum extrapolations are 
X 2 /Ndf ~ 0.5 for mesons, « 1.0 for octet baryons and « 
3.0 for decuplet baryons. 



B. QxPT fits in the continuum limit 



We fit the continuum hadron spectrum, normalized by 



(0.75) 
V 



using the QyPT formulae and following the pro- 



Tfiy , using the U.yr J 

cedure given in Sec. [Vj, and obtain the hadron spectrum 
as well as the value of m { y 75) = 0.981(18) GeV. The 
quark mass in the continuum limit at each mps/my is 
necessary for a chiral fit of pseudoscalar meson masses. 



AWI(0) 



lin- 



We first interpolate the AWI quark mass 
early in 1/k to each value of mps/rny at each (3. We 
then convert rriq WI ^ to the values in the MS scheme 
at n = 2 GeV, using a determined from mj, 0,75 '. The 
AWI quark masses are linearly extrapolated to the con- 
tinuum limit with a reasonably s mall x 2 /Ndf < 0.5. For 
completeness, we provide in Table XL! lattice spacings at 



each f3 and quark masses for each mps /my normalized 

, (0.75) 

by m v . 

The QxPT fits in the continuum limit look quite simi- 
lar to those at finite a as shown in Fig. The values of 
the parameters are consistent with those obtained from 
an extra polat ion of the parameters at finite (3, as given 
in Tabic XLII . See also, for example, Fig. ^5] in which we 
plot the octet baryon mass in the chiral limit. A compar- 
ison of the coefficients of the singular terms of the QxPT 
formulae is made in Appendix ^. 



15 



C. Universality of the light hadron spectrum 



Table XLID summarizes the light hadron masses thus 



obtained together with deviations from those from the 
original procedure. The two spectra are consistent with 
each other with differences smaller than la of the statis- 
tical error. See also Fig. ^ in which we compare masses 
from the two methods for the case of the uik input. For 
interpolation/extrapolation of hadron masses, we test 
polynomial chiral fits instead of the QxPT formulae. The 
deviations remain within la for pseudoscalar mesons and 
baryons and 1.5a for vector mesons. 

We therefore conclude that the two spectra determined 
by the two methods are consistent with each other with 
differences smaller than 1.5a. This confirms that both 
chiral and continuum extrapolations are under control. 
We take the differences as systematic errors due to the 
chiral and continuum extrapolations. 



XI. CONCLUSIONS 

In this article we presented details of our calculation of 
the light hadron spectrum and quark masses in quenched 
QCD. The computational power provided by the CP- 
PACS computer enabled an exploration of hadron masses 
at lighter quark masses than hitherto attempted. 

The high precision data for pseudoscalar meson masses 
revealed evidence supporting the presence of chiral sin- 
gularities as predicted by quenched chiral perturbation 
theory. In the vector meson and baryon sectors the preci- 
sion of our data is not sufficient to draw conclusive state- 
ments on quenched chiral singularities. However, simula- 
tions covered the range of quark masses sufficiently small 
to obtain a stable result of the spectrum also for vector 
mesons and baryons. Predictions do not depend on the 
choice of conventional polynomial chiral fits or fits based 
on quenched chiral perturbation theory. Since scaling 
violation also turned out to be mild for the plaquette 
gluon and Wilson quark actions, the hadron masses keep 
the statistical precision of 1-2% for mesons and 2-3% for 
baryons after the continuum extrapolation. The system- 
atic error is estimated to be at most 1.7%. 

The chief finding is the pattern and magnitude of the 
breakdown of the quenched approximation for the light 
hadron spectrum. In the meson sector the quenching 
error manifests itself in a small hyperfine splitting when 
compared with experiment. A small mass splitting is 
also seen in the decuplet baryons, and masses themselves 
are small for octet baryons. The magnitude of deviation, 
typically 5-10%, is much larger than the statistical and 
systematic errors. 

The quenched approximation poses a limitation of our 
ability to predict fundamental parameters of QCD. The 
strange quark mass m s depends on hadron mass input, 
with a difference as large as 25%. The QCD scale pa- 
rameter has uncertainty of the order of 15% depending 



on inputs, i.e., m p or phenomenological values of y/a or 
ro- 

It appears to us that it is not worthwhile to pursue 
precision further, and the effort should rather be directed 
toward QCD simulations incorporating the effects of dy- 
namical sea quarks. We in fact started an attempt |5l| 
in this direction using improved gluon and quark actions. 
In the course of this work we recalculated the quenched 
hadron spectrum using the improved actions. The light 
hadron spectrum obtained in the continuum limit is in 
good agreement with the results reported in this article, 
providing further confirmation of the success and limita- 
tion of quenched QCD. 
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APPENDIX A: GAUGE FIXING 

For the measurement of hadronic observable, we fix 
gauge configurations to the Coulomb gauge by maximiz- 
ing the quantity 



(Al) 



11=1 



To achieve this, we combine two methods: (a) an SU(2) 
subgroup method, which is similar to the pseudo heat- 
bath algorithm for SU(3) gauge theories pl|, and (b) an 
over-relaxed steepest descent method |32| . Both methods 
can be vectorized and parallelized by splitting the lattice 
sites into even and odd sites. 



1. SU(2) subgroup method 



Under a gauge transformation 



(A2) 



with SU(3) matrices G n , H transforms as 
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H -> H' = 1* [GnC/,^ + 17, 
+terms indep. of G 



(A3) 



If the gauge group is SU(2), it is easy to find the solution 
G n which maximizes H' for a given site n. The global 
maximum can be achieved iteratively by repeating the 
maximization at all sites. For the SU(3) case, we can 
gradually increase H by applying the maximization for 
different SU(2) subgroups of SU(3). In our simulations, 
we maximize three SU(2) subgroups per iteration. 



configuration is far from the maximum, however, this 
method sometimes fails to converge. Therefore, we first 
apply the SU(2) subgroup method for several hundred 
iterations to drive the configuration close to the maxi- 
mum. In our simulations, we adopt the SU(2) subgroup 
method for the first 200, 500, 1000 and 6000 iterations at 
/3 =5.9, 6.1, 6.25 and 6.47, respectively, before applying 
the steepest descent method. 

For a convergence check, we monitor 



h = H/(9L 3 s L t ) 



(All) 



and 



2. Over-relaxed steepest descent method 

Alternatively, H can be maximized by iterative gauge 
transformations with 



G n = exp(a„A„) 

3 

A n — ^ ^ A n; ^ 
u=l 



- U, hfJl - h.c. - trace] 



(A4) 
(A5) 
(A6) 



with a suitable real number a n , because the gradient of 
H with respect to 6® defined by G n — exp(iA a 6*°) is given 
by 



SH 

— = ReTr[-iA°A r , 



(A7) 



Note that vanishes at the supremum of H and 

is proportional to £ d^A^ix) in the continuum limit. 

A candidate of a n is obtained by solving the supremum 
condition of H 



Re^ Tr[e Q " A "A„C/„, At - C/„_ AiM A„e- Q " A "] = (A8) 

u=l 

to the leading order of a n , which gives 

_ R e E M Tr [A n {C/ n _ A)/1 - U n ^}] 



Re£„Tr [& n {U n ^ + U n , M }] 



(A9) 



An over-relaxation is introduced through a parameter 
1 < lu < 2 in the gauge transformation: 

G n = exp(wa„A„). (A10) 



3. Implementation on the CP-PACS 

We find that the steepest descent algorithm with the 
over-relaxing parameter to w 1.98-1.99 converges much 
faster than the SU(2) subgroup method, when the con- 
figuration is already close to the maximum. When the 



^ = i^E ilr[A n , M Atj. (A12) 



Note that 



L 3 L t ^3^3 

s 1 n /i=l 



. 3 8 

Acx / dxJ2J2( d » A l 

J n=l a=l 



(A13) 



in the continuum limit. We truncate iterations at z-th 
iteration when the conditions 



\hi -hi-!] < 10- 10 , 
A, < 10- 14 , 



(A14) 
(A15) 



are both satisfied. We have checked that stronger con- 
vergence criterion does not lead to a significant difference 
in hadron propagators. 



APPENDIX B: QUARK PROPAGATORS 



In order to solve Eq.(17), we use a red/black- 
preconditioned minimal residual (MR) algorithm. We 
accelerate the convergence by applying successive over- 
relaxations. For the over-relaxation factor, we adopt 
/ = 1.1 from a test study made at / = 0.9 1.0, 1.1 
and 1.2. The number of iterations for / = 1.1 is smaller 
than those for / = 0.9 (/ = 1.2) by 20% (5%). 

At each MR step, we monitor the residual sum of 
squares R 2 , where B = S — D(k)G, and truncate iter- 
ations when the condition R 2 < 10~ 12 (R 2 < 10~ 14 at 
j3 = 6.47) is satisfied for the point source, and R 2 < 10~ 7 
for the smeared source. Hadron propagators obtained 
with this stopping condition are compared with those 
with much stronger one on several configurations. From 
this test we estimate that the truncation error in hadron 
propagators on each configuration is smaller than 5% of 
our final statistical error for any particle at any time slice. 

The number of iterations needed to calculate quark 
propagators is listed in Table XLIV . The number is ap- 
proximately proportional to the inverse of the quark mass 
defined by 1/k c )/2, where k c is the critical hopping 

parameter. 
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Our exponentially smeared source Eq. fllq ) is moti- 
vated from a result of the JLQCD Collaboration for the 
pion wave function, 



*(r) 



^T.n^i^^jn + r)\Tr) 



(Bl) 



which was well reproduced by a single exponential func- 
tion \I>(r) = Acxp(— Br) except at the origin ^(O) w 1.0 
plf. The coefficient A and the slope B of the JLQCD 
Collaboration can be parametrized as 

A(m p a, m q a) — a,Q + a\m p a + (02 + a 3 (m p a) 2 )m q a, (B2) 
B(m p a, m q a) — b + bim p a + (b 2 + b 3 (m p a) 2 )m q a, (B3) 

where m p a is the dimensionless p meson mass in the chi- 
ral limit, m q a = (1/k — 1/k c )/2 is the bare quark mass, 
and 

a = 0.915, ai = 0.576, 
a 2 = 0.2127, a 3 = -0.644, 
b = -0.0537, 61 = 0.978, 
b 2 = 0.2146, 63 = -0.5123. 

Applying the results of test runs for m p and k c , we adopt 
A and B listed in Table KLIV. The smearing radius is 
approximately constant, a/B 0.33 fm at our four j3 
values. 



APPENDIX C: HADRON PROPAGATORS 

For mesons we employ the operators defined by 

M f A 9 (n)=J n T A g n (CI) 

where / and g are quark fields with flavors / and g, and 
r^4 is one of the 16 spin matrices 

T s = I, T P = 75, Tp = ^7075, IV = H, I> = iloli 
T A = ij 5 ^ and T T = i[ 7i , 7j ]/2 (i, j = 1, 2, 3). (C2) 

With these operators, we calculate 16 meson propagators, 
(M A (n)M A (0)). 

For the spin 1/2 octet baryons we take the operators 
defined by 



0£"» = e abc (f^C l5 g b n )h c nl 



(C3) 



where a, 6, c are color indices, C = 7472 is the charge con- 
jugation matrix, and a = 1,2 represents the spin state, 
up or down, of the octet baryon. To distinguish S- and 
A-like octet baryons, we antisymmetrize flavor indices, 
written symbolically as 



E = — 
A = 



[fh]g+[gh]f 

V2 ' 
[fh]g-\gh]f-2[fg]h 

^/6 



(C4) 
(C5) 



with [/s] = 1 9 - gf- 

The spin 3/2 decuplet baryon operators are given by 
Dl 9 » = t abc (fn a C^9 b n )h c na . (C6) 
Writing out the spin structure (fx, a) explicitly, we obtain 

C3/2 = e a " c (/ Ta CT +5 b )^, (C7) 

D 1/2 = e abc [(f Ta CT og b )hl - (f Ta CT+g b )h c 2 }/3, (C8) 

D-i/2 - e abc [(f Ta Cr g b )h c 2 - (/ Ta CT_.g h )^]/3, (C9) 

D^ 3/2 = e abc (f Ta Cr_g b )h c 2l (C10) 

where T± = (71 =F i~f 2 )/2, Tq — 73, and the subscript of 
D denotes the z-component of the spin. 

With these operators, we calculate 8 baryon propaga- 
tors given by 

(S Q (n)S Q (0)), a = 1,2 (Cll) 
(A Q (n)A a (0)), a =1,2 (C12) 
(D s (n)D 8 (0)), S = 3/2,1/2,-1/2,-3/2 (C13) 

together with 8 antibaryon propagators defined by the 
same expressions with the baryon operators replaced by 
antibaryon operators. To enhance the signal, we aver- 
age zero momentum propagators on a configuration over 
all states with the same quantum numbers; three polar- 
ization states for the vector meson and two (four) spin 
states for the octet (decuplet) baryon. We also average 
the propagators for the particle and the antiparticle, i.e., 
meson propagators at t and L t — t are averaged, and 
baryon propagators for particle at t and those for an- 
tiparticle at L t — t are averaged. Errors of propagators 
are estimated treating the data thus obtained being sta- 
tistically independent. 



APPENDIX D: CORRELATED FITS FOR 
CHIRAL EXTRAPOLATION 

A difficulty in a correlated chiral extrapolation is that 
the size of the full covariance matrix (error matrix) 
C(t, K\K 2 \ t' , n'iK 2 ) is very large and the matrix becomes 
close to a singular matrix so that C~ 1 necessary for 
X 2 fits cannot be estimated reliably. When we make 
a fit for both degenerate and non-degenerate data si- 
multaneously, the size of C becomes of order 200, e.g., 
(28 - 10 + 1) x 11 = 209 for fitting range [10,28] used 
for 11 combinations of quark masses. We find that the 
condition number of C is far beyond 10 15 so that one 
cannot handle the matrix within our numerical accuracy. 
Instead of the simultaneous fits, we make independent 
fits for degenerate and non-degenerate cases. Namely we 
make three fits for mesons: 1) degenerate fit with 5 data 
at K\K 2 = s±si, s 2 s 2 , U1U1, u 2 u 2 , U3U3, 2) non-degenerate 
fit with one of the hopping parameter k\ — s\ using 4 
data at k,ik 2 = sisi,s\ui,siu 2 ,siu 3 , and 3) the same 
with Ki = s 2 using kik 2 = s 2 s 2 , s 2 U\, s 2 u 2 , s 2 u 3 . Fits for 
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baryons are carried out similarly, because two quarks in 
baryons are taken to be degenerate. 

For each correlated fit, we employ the procedure 
adopted in Ref. |53| . We first minimize x%u defined by 



Xfull 



{G(t, k 2 ) - G (t, n 2 )}C-\t, K2 ;t' K ' 2 ) 

{G(t', K ' 2 )-G (t', K ' 2 )}, (Dl) 



where G(t, k 2 ) are data of hadron propagators and 
Go(t, k 2 ) is a fitting function e.g. Go(t,n 2 ) = 
A(k 2 ) exp(— m(K 2 )t) for baryons. The masses tti{k 2 ) thus 
determined arc in general different from those obtained 
by individual x 2 fits f° r each n 2 . The difference is small 
for most cases, though it occasionally amounts to 1.2a. 
We use masses from the full correlated fits for later anal- 
yses. Results remain essentially the same if masses from 
individual fits are used. We then calculate an error ma- 
trix £ for the fit parameters by 



S = (D T C- 1 D)- 1 , 
where D is the Jacobian defined by 



(D2) 



= [dG (t, K 2 )/dA( K ' 2 ),dG (t, « 2 )/9m(4)]. (D3) 

Note that D is diagonal with respect to n 2 ■ For the chiral 
extrapolation, we minimize x\xt given by 



for the former, we observe that X%u/Ndf for the latter are 
much larger than unity although errors of X%u/Ndf are 
also large. After trial and error, discarding data around 
t w t m i n in degenerate propagators at K = Si and s 2 
and/or u\ leads to X%u/Ndf ~ 1- We therefore use differ- 
ent ranges for polynomial chiral fits from those for QxPT 
fits, noting that masses for these cases do not change sig- 
nificantly. 

For the mass-mass covariance matrix E, the condition 
number is O(10 4 ) and all errors for eigenvalues and eigen- 
vectors are contained within 25% of central values. Hence 
we are able to perform numerically reliable full correlated 
chiral extrapolations. 



APPENDIX E: CHIRAL FITS FOR 
PSEUDOSCALAR MESONS 

We compare eight chiral fi t func tions for pseudoscalar 
meson masses listed in Tabic XLV1 , using fully-correlated 
fits described in Appendix [5|. The first five are for the 
degenerate cases, SiSi and UiUi. The fits 1-3 are polyno- 
mials, while the fits 4 and 5 are based on QxPT using 
Eq. (|2^) with a$ = 0, with or without the quadratic 
term B(mi + m 2 ) 2 . The remaining three functions are 
for the non-degenerate cases, SiUj. Two of them are poly- 
nomials, and the last one is the QxPT formula ( f28| ) with 
ctq, = 0. For non-degenerate fits, we fix n c to a value 
determined from a degenerate fit. 



Values of Xext/^df for chiral extrapolations are very 



X-ext 



= y^{m{K 2 ) - /(K 2 )}E _1 (K 2 ,K2){ m ( K 2) - f( K ' 2 )}> large irrespectiv e of t he choice of fitting functions 



(D4) 

where f{n 2 ) is the fitting function we try and the matrix 
T.(k 2i k 2 ) is the sub-m atri x among the masses of the full 
error matrix £ in Eq.( |D2| ). 

This procedure works well only when the full covari- 
ance matrix C is reliably determined. Although the con- 
dition number of C is as large as 10 14 - 10 15 for degenerate 
fits and 10 12 -10 13 for non-degenerate fits, small eigen- 
values and the corresponding eigenvectors responsible for 
Xfuii are determined well. The jackknife error for these 
quantities is of the order of 10% (20% at f3 = 6.47). The 
error sometimes increases to 50% for large eigenvalues. 
We think that it causes no problem because correspond- 
ing eigenmodes have no significant contribution to X%u- 

In fully correlated chiral fits, x%a should be close 
to Ndf- For pseudoscalar mesons, we choose the fit- 
ting range carefully to satisfy this condition (see Table 



shown in Table XLVI . A similar phenomenon was ob- 
served also in previous studies. See e.g. Ref. p). This 
may be due to the fact that higher order terms are re- 
quired to reproduce our data. Because the number of our 
data points is limited, inclusion of such terms is not pos- 
sible, and hence we choose a functional form from overall 
consistency. 

Concerning the relative magnitude of Xext/^df, we 
find, for the degenerate fits, that Xext/^df is the small- 
est with the QxPT formulae keeping the 0(m 2 ) term 
(Fit 5). When we remove the 0{m 2 q ) term, xlxt/^df be- 
comes much larger. This observation is consistent with 
the presence of the (3(m?) term expected from the mass 
ratio test given in Sec. M. For non-degenerate fits, simi- 
lar values of x\xtl^df are obtained from both quadratic 
(Fit 2n) and QxPT (Fit 5n) fits. 

To keep consistency with the presen ce of QxPT singu- 
larity shown by the ratio tests in Sec. VA, we decide to 



XLV). We use the common ranges for both Qx PT fit s in 
Sec. VIA and quadratic polynomial fits in Sec. VIE. 

For vector mesons and baryons, we employ uncorre- 
cted QxPT chiral extrapolations made simultaneously 
to hadron masses with degenerate and non-degenerate 
quark masses. In addition, we perform fully correlated 
polynomial chiral fits. With our choice of fitting ranges 



employ QxPT fits (Fit 5 and Fit 5n) for the main course 
of our analyses and use quadratic fits (Fit 2 and Fit 2n) 
for estimations of systematic errors from the chiral ex- 
trapolation. 
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APPENDIX F: TEST OF QxPT MASS FORMULA 
FOR VECTOR MESONS AND BARYONS 

Lowest order QyPT mass formulae for vector mesons 
jl7j and baryons |lq] can be written as 

m H (mp S ) = m° + C 1/2 mp S + Cim 2 PS + C s / 2 m PS , 

(Fl) 

where Cj are polynomials of the couplings in the 
quenched chiral Lagrangian. We find that it is difficult 
to constrain all coupling parameters (6 for vector mesons 
and 11 for baryons in addition to S and a$) under the 
limitation of the accuracy of our mass data and the num- 
ber of data points. We, therefore, set S — 0.1 and a$ = 0, 
and drop the couplings of the flavor-singlet pseudoscalar 
meson to vector mesons and baryons. We also set / = 132 
MeV unless otherwise stated. 



1. Vector Mesons 



U = Ci/2 + C lX + 0{m 2 PS ) 



for 



y 



my,22 - m v ,ii 



mps,22 - rnps,n 
x = mps,22 + mps,u- 



(F7) 

(F8) 
(F9) 



We c alcu late y and x for all 10 combinations of m and k 2 . 
Eq. (F7) is obtained also for non-degenerate cases, with 
y and x replaced by more complicated expressions. We 
obtain 15 data points for y and x from all combinations 
satisfying (uiSj) ^ (m's' ). 

In Figs. [|6] and |47] we show plots of y versus x for 
degenerate and non-degenerate cases, respectively. Data 
for y are fitted well by a linear function of x and intercepts 
are negative taking a value in the range —0.3 - 0.0. These 
results suggest that the 0(mp s ) term in Eq . ( |F7|) and 
hence C 3 / 2 in Eq.(|F6|) or Cd and Cat in Eq.(|F2j) are 
small. We find that C\/2 is negative but much smaller in 
magnitude than the phenomcnological estimate ss —0.71. 



The lowest order Q%PT formula for vector mesons |17 
is given by 



b. Chiral fit 



C 3 i 
m v = m v H ]^-{^(m uu + m ss ) + 2- 

C\ . 9 2 \ 

+ -^-("V U + m ss) 

+C D (ml u + m 3 ss ) + C N m 3 us , 



-} 



(F2) 



where m/ s is the pseudoscalar meson mass. The coef- 
ficients are written in terms of the couplings gi of the 
quenched chiral Lagrangian; 

Ci /2 = -47r 32 2 «5, (F3) 
C D = -2 3254 /(12^/ 2 ), (F4) 
C N = (-4g l92 + 4g 2 2 A N )/(127rf). (F5) 

The coefficient Gyn of the term linear in mps is pro- 
portional to S, and represents the quenched singular- 
ity- QxPT predicts a negative value for Ci/ 2 . A phe- 
nomenological estimate is Ci/2 ~ —0.71 using S = 0.1 
and g 2 = 0.75 pi 



a. Ratio test 

We perform ratio tests for vector meson masses inde- 
pendently for degenerate and non-degenerate cases. For 
the degenerate case, the mass formula Eq. ( |F2] ) reduces 
to Eq. (plh with 



C3/2 — Cu + Cat. 
Hence, we obtain a relation 



(F6) 



We make a fit (F2) directly to the vector meson mass 
data, treating the degenerate and non-degenerate cases 
simultaneously. We ignore correlations among masses for 
different quark masses, or else the size of full covariance 
matrix becomes too large to obtain reliable matrix ele- 
ments. 

We find that the QxPT fit keeping all five fitting pa- 
rameters is unstable; the covariance matrix for the fit 
parameters becomes close to singular with the condition 
number of O(10 8 -10 9 ). 

Dropping the 0(m PS ) terms, the fit become more sta- 
ble with the condition numbers of C(10 4 -10 6 ). The fit re- 
produces data equally well as that including the 0(mp s ) 
terms as illustrated in Fig. ^7| for the degenerate data 
at j3 = 5.90 (the baryon fits in this figure are discussed 
below). Equivalently, x 2 /Ndf of at most 0.8 obtained 
without the 0(mp s ) terms are comparable to 0.9 includ- 
ing the 0(mp s ) terms. Taking the stability of fits as 
a guide, we adopt the fit without the 0(m PS ) terms for 
vector mesons. This choice also agrees with a small value 
of C3/2 observed in the ratio test. 

For the coefficient of the leading chiral singularity, we 
obtain C 1/2 = -0.058(27), -0.075(28), -0.065(35) and 
-0.155(72) at = 5.90, 6.10, 6.25 and 6.47, respectively, 
which are stable under variation of (3. Taking a weighted 
average, we find C1/2 = —0.071(8) which is much smaller 
than a phenomenological estimate C1/2 = —Aitg^S — 
-0.71 (5 = 0.1). The value C 1/2 = -0.13(10) obtained 
from chiral fits in the continuum limit (Sec.|x]) also sup- 
ports this conclusion. A much larger value of C1/2 ~ 

—0.6 0.7 is reported in Ref. |^9|. One possible origin 

of the difference is finite size effects in Ref. If29f . 
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2. Decuplet baryons 

Keeping terms up to 0(mp s ), the QxPT formula for 
decuplet baryon masses is given by Q 

5H 2 

m D = m° D + ——{Awuu + 4w ua + w ss ) 
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C 2 



+ 1~8 ^ Wuu ~ 2Wus + Wss ^ 

+ c ( 2m lu + m ls) 

+ (-10ff 2 /81 + C 2 /9){v uu + 2v us ), (F10) 

where Wij and Vij are defined in Eqs. ( [43| ) and (^) , The 
terms proportional to Ujj are 0(mp s ). 



a. Ratio test 



For the degenerate case, the formula (F1C) reduces to 
the cubic polynomial as in Eq.(Fl), where 



Ci/2 = -(5irS/6)H 2 



(Fll) 



Therefore, we can perform a ratio test by constructing y 
and a: which are similar to those for vector mesons. 

As Fig. |^| shows, we observe C\i-x ~ —0.1 from the 
intercept and a small value of C3/2 from the linearity of 
data. The negative value for Ci/ 2 is consistent with the 
negative sign in Eq.(Fll). 



b. Chiral fit 

We study direct fit s of decuplet baryon mass data 
to the mass formula (|F10| ) using degenerate and non- 
degenerate masses simultaneously. Since C3/2 is small 
in the ratio test, we consider fits with and without the 
0(m PS ) terms in Eq. flFlO] ). We note that the number of 
parameters is the same (four). 

The two types of fits are indistinguishable. Both yield 
X 2 /Ndf ~ 0.5 with the condition number of the covari- 
ance matrix of O(10 6 -10 7 ), and the fitting curves are 
nearly identical in the range of our data points. See 
Fig.|27] for the results at /3 = 5.90. 

The fits without the 0(mp s ) terms give H 2 = 
0.65(13), 0.55(15), 0.39(15) and 0.70(40), and C 1/2 = 
-(5nS/6)H 2 = -0.169(33), -0.145(39), -0.101(40) and 
-0.18(10) for f3 = 5.90, 6.10, 6.25, and 6.47, respectively. 
Results at different (3 are consistent with each other, and 
the weighted average C]/ 2 = —0.14(1) is consistent with 
the estimate sw —0.1 from the ratio test. 

Including the 0(m PS ) terms reduces the Ci/ 2 coeffi- 
cient to C 1/2 = -0.020(29), -0.019(32) -0.031(28) and 
-0.204(73) for (3 = 5.90, 6.10, 6.25, and 6.47, which are 
much smaller than the values from the ratio test. 

Considering the consistency with the ratio test, we em- 
ploy the mass formula without the 0(m PS ) terms in the 



main analyses, and use the fits with the 0(m PS ) terms 
for error estimations. 

For the adopted fit, C 2 is consistent with zero within 
~ 2a. Setting C 2 = does not change the decuplet 
baryon masses at the physical point by more than 0.5a 
for all f3 values. We keep, however, the term proportional 
to C 2 since it is the same order as the term proportional 
to H 2 . 



3. Octet baryon 

The QxPT mass formulae we consider for octet 
baryons are written as 



ms = m 

+ ^{4F 2 w uu - A{D - F)Fw us + (D - F) 2 w ss } 
C 2 

+-^-{w uu - 2w us + w ss ) 

-ib F m 2 uu + 2(b D - b F )m 2 ss 
+ (2D 2 /3-2F 2 -C 2 /9)v uu 
+ {2D 2 /3-4DF + 2F 2 -5C 2 /9)v us , (F12) 
toa = rn o 

+ i{(4TJ/3 - 2F) 2 w uu + (D/3 + F) 2 w ss 

-2(4D/3-2F){D/3 + F)w us } 

+4(26^/3 - b F )m 2 uu - 2{b D /Z + b F )m 2 ss 

+ {2D 2 /9 - 8DF/3 + 2F 2 - C 2 /3)v uu 

+ (10L> 2 /9 - ADF/i - 2F 2 - C 2 /3)v us . (F13) 

The notation is the same as that for decuplet baryons. 



a. Ratio test 



The QxPT mass formula (F12) reduces to the cubic 
polynomial in Eq. (Fl) with 



Cl/2 



-(3tt<5/2)(L>-3F) :; 



(F14) 



The negative sign of Gxn suggests that the degenerate 
octet mass is a concave function of mp s for sufficiently 
light quarks. On the other hand, our data exhibits a 
convex (i.e., negative) curvature as shown in Fig. A 
negative curvature is also observed in Ref. || for the 
Kogut-Susskind action. We consider that the negative 
curvature is due to 0(mp s ) terms, and that they have a 
large effect for the range of quark masses covered by our 
data. 

In Fig. [49 we show the ratio test as in Eq. (F7) for de- 



generate quark masses. If 0(m 2 PS ) terms in Eq. (F7), or 
equivalently the 0(mp s ) terms in the mass formula, are 
negligible, one would obt ain C 1/2 ~ 0.6 which is opposite 
in sign compared to Eq. (F14). 
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b. Chiral fit 

For octet baryons, it is natural to fit S- and A-likc 
baryons simultaneously, because many coefficients of in- 
dividual terms are related with each other. In these fits 
we need to include the 0(m PS ) terms to reproduce the 
negative curvature while maintaining consistency with 

QxPT. 

We find that the six-parameter fit as indicated by 
Eqs. (F12) and (F13) shows several local minima in pa- 
rameter space, and the minimization procedure converges 
to different minima depending on the jackknife ensemble. 
Some parameters exhibit irregular dependence on the lat- 
tice spacing, e.g., F = 0.150(20), 0.139(25), 0.103(43) 
and 0.247(49) for /? = 5.90, 6.10, 6.25, and 6.47. 

Recalling that the octet-decuplet coupling C is zero 
within 2a in the decuplet baryon fit, we set C = for 
octet baryons and obtain stable fits, e.g., F = 0.334(14), 
0.326(14), 0.315(13), 0.299(31) for (3 = 5.90, 6.10, 6.25 
and 6.47. 

In view of the stability of parameters, we employ the 
simplified fit (C = 0). The fits are almost indistinguish- 
able from those keeping C ^ in the range of measured 
points. See Fig.p7[ 

For the coefficient of the leading chiral singularity, we 
obtain Cj/2 = —0.118(4) as the weighted average over the 
four values of f3. This is smaller than a phenomenological 
estimate, C 1/2 = -(3nS/2)(D - 3F) 2 w -0.27 assuming 
S = 0.1, F = 6.5 and D = 0.75. The smallness of C 1/2 is 
not due to lattice cut-off effects, because C±n = —0.09(8) 
is obtained in the continuum limit (Sec.[x|). 



APPENDIX G: QUARK MASS FROM LOCAL 
AXIAL- VECTOR CURRENT 

An alternative definition of the AWI quark mass is 
given by 



where tops a is already determined by the pseudoscalar 



2m £ WI > local a 



ylocal 

^_ 2 AWI,loca m ^ (G1) 

Zp 

O^AWI,local(0) _ „ i- {A l £ Cal (t)P(0)) 

" ^^Too (P(i)P(O)) ' (G2) 

where a time derivative in the numerator is substituted 
by tops a, and 



4 ocai (n) = /„i 7 57^ 
is the local axial- vector current. 



(G3) 



For to 



AW I ,local(0) 



we first antisymmetrize the correla- 



tor (A^ cai {t)P(0)) and then make a fit 



rnpsa- 



(Af cal (t)P(0)) 



<P(i)P(0)) 

2m Awi,local(Q) tanh(TOpsa(L t /2 - t)), (G4) 



propagator fit . The r esults for 2m 



AWIlo 



22(0) 



rized in Table XLVII 



polynomial under the constraint 

the renormalization coefficient for m 
Eq. p|). 



are summa- 
Chiral fits are made by a quadratic 

AWI,local _ ^ p or 

AWI,local(0) 
q ' ) 



we use 
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TABLE I. Simulation parameters. The lattice spacing a is determined from m p . The last column "acceptance" is the mean 
acceptance rate in the pseudo-heat-bath update sweeps. 



p 


Li x L t 


o-^GeV] 


a[fm] 


L s a[fm] 


#conf. iter./conf. 


acceptance 


5.90 


32 3 56 


1.934(16) 


0.1020(8) 


3.26(3) 


800 200 


0.85 


6.10 


40 3 70 


2.540(22) 


0.0777(7) 


3.10(3) 


600 400 


0.84 


6.25 


48 3 84 


3.071(34) 


0.0642(7) 


3.08(3) 


420 1000 


0.83 


6.47 


64 3 112 


3.961(79) 


0.0498(10) 


3.18(6) 


150 2000 


0.82 






TABLE II. 


Hopping parameters and corresponding 


mps/mv ratios. 












P = 5.90 






K 




1 5660 


0.15740 


0.15830 


1 5890 


1 5920 


, 1 1 . , ^, /mo t r 
'HPS 1 IfiV 




yj. I tJZi ^ L j 


0.692(1) 


0.593(1) 


KJ.'-tiJ -L I ^ J 












/3 = 6.10 






K 




0.15280 


0.15340 


0.15400 


0.15440 


0.15460 


mps/mv 




0.751(1) 


0.684(1) 


0.581(2) 


0.474(2) 


0.394(3) 










/3 = 6.25 






K 




0.15075 


0.15115 


0.15165 


0.15200 


0.15220 


mps/mv 




0.760(1) 


0.707(2) 


0.609(2) 


0.502(2) 


0.411(3) 










P = 6.47 










0.14855 


0.14885 


0.14925 


0.14945 


0.14960 


mps/mv 




0.760(2) 


0.709(3) 


0.584(3) 


0.493(4) 


0.391(4) 






TABLE III. Measured execution time of each step 


in units of hour. 




P 






5.90 


6.10 


6.25 


6.47 


size 






32 3 x 56 


40 3 x 70 


48 3 x 84 


64 3 x 112 


#PU 






256 


512 


1024 


2048 


configuration generation 




0.28(10%) 


0.65(14%) 


1.62(24%) 


4.5(29%) 


gauge fixing 






0.13(5%) 


0.22( 5%) 


0.33(5%) 


1.8(12%) 


quark propagator 




1.57(52%) 


2.72(57%) 


3.13(46%) 


6.6(42%) 


hadron propagator 




1.01(33%) 


1.18(24%) 


1.67(24%) 


2.6(17%) 


total 






2.99 


4.76 


6.78 


15.6 
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TABLE IV. Quenched light hadron mass spectrum. The first error is statistical and the second is systematic. Deviation 
from experiment and its statistical significance are also given. 



liei axon 


6xp GrimGirt 


ttlk input 




uifi input 






TP r>Vl 


mass [Ge V] 


deviation 


mass i Vjre v I 


deviation 


t\ 


A /1Q77 






u.ooo^ iu ) y\Jo ) 


1 1 . z /o o . Da 


jy'* 
K 


u.oyoi 


0.808(09^08) 


— 4.27o 4.ja 


0.889(03) (06J 


—0.87o 2.3a 


<p 




O OK^I /'I A\ 

u.yo I (Ld) (i-'i) 


(? 10? A Q _ 

— O.l/o 4. OCT 






N 


0.9396 


0.878(25) (14) 


-6.6% 2.5a 


0.878(25)(14) 


-6.6% 2.5a 


A 


1.1157 


1.019(20) (09) 


-8.6% 4.7cr 


1.060(13) (10) 


-5.0% 4.1a 


£ 


1.1926 


1.117(19)(11) 


-6.4% 4.1a 


1.176(11)(20) 


-1.4% 1.5a 




1.3149 


1.201(17)(13) 


-8.7% 6.8a 


1.288(08) (09) 


-2.0% 3.5a 


A 


1.2320 


1.257(35)(10) 


2.0% 0.7a 


1.257(35) (10) 


2.0% 0.7a 


S* 


1.3837 


1.359(29)(11) 


-1.8% 0.9a 


1.388(24)(11) 


0.3% 0.2a 




1.5318 


1.459(26)(10) 


-4.7% 2.8a 


1.517(16)(09) 


-1.0% 0.9a 


Q 


1.6725 


1.561(24) (09) 


-6.7% 4.7a 


1.647(10)(15) 


-1.5% 2.6a 



TABLE V. Decay constants for light pseudoscalar and vector mesons. 







experiment 


quenched QCD 


U [MeV] 




130.7(0.4) 


120.0(5.7) 


f K [MeV] 


(m-K-input) 


159.8(1.5) 


138.8(4.4) 




(ra^-input) 




141.5(3.8) 


U/Ik - 1 


(mx-input) 


0.223(2) 


0.156(29) 


F p [MeV] 




220(5) 


205.7(6.6) 


F<t, [MeV] 


(m^-input) 


239(3) 


229.4(5.7) 



25 



TABLE VI. Meson mass data in lattice units at the simulation points. The fit range [tmin,tmax] and X 1 /Ndf are also listed. 
Errors of \ /Ndf are estimated by the jackknife method. 



P = 5.90 







Pseudo- Scalar s 






Vectors 






range 


X /Ndf 


mass 


range 


2 / AT 

X /Ndf 


mass 




10 28 


0.55(36) 


0.38255(25) 


7 28 


1.03(46) 


0.50900( 60) 


S2S2 


10 28 


0.49(34) 


0.33114(26) 


7 26 


1.11(51) 


0.47862( 71) 


Ui til 


1 n ofi 
1U zo 


U.4y^o4 ) 


U.ZD41 1 ^Zo J 


7 OO 

/ zz 


U.o4^oz ) 


U. 44014^ yDJ 


U2U2 


1 n 98 


U.OO^OD ) 


0.20827(33) 


7 18 


U.oO^oo ) 


fl 49^01 ( 1 39"! 

u.4zoyi^ioz j 




1 n os 

1U Zo 


U. < 1 yoo ) 


0.17145(54) 


7 16 


C\QAf QQ\ 

u.o4^oy ) 


u. 41011^1 / OJ 


<5 1 u 'l 


10 28 


50(34^ 


0.32833(26) 


7 26 


1.11(51) 


47749 ( 741 


SlU2 


10 28 


51 (351 

W . U 1 1 tJ'J 1 


0.30769(28) 


7 24 


1 (Wfv^ 

-L . VJ ij I <J tj J 


46724( 841 


1 Lt u> 


10 28 


51(351 


0.29666(31) 


7 22 


34(321 


46230( 931 


S2U1 


10 28 


49(341 


0.29945(27) 


7 24 


1 08(541 

-L . VJU I <Jt: 1 


461 98( 821 


S2U2 


10 28 


49(341 


0.27674(28) 


7 22 


30(301 


451 62( Q31 




10 28 


0.47(34) 


0.26440(31) 


7 21 


35 f 341 


44666(1031 








(3 = 6.10 










Pseudo-Scalars 






Vectors 






range 


X /Ndf 


mass 


range 


2 / AT 

X /Ndf 


mass 


SlSl 


14 35 


1.49(57) 


0.29143(22) 


9 35 


0.76(36) 


0.38788( 53) 


S2S2 


14 35 


1.51(58) 


0.24706(23) 


9 32 


0.57(33) 


0.36144( 64) 


tiitii 


14 00 


1 ^n/^7^ 
l.OU^O ( ) 


0.19514(24) 


9 26 


f\ on/ oo\ 


fl QQCQQ/ QO\ 


U2U2 


14"^ 




0.15170(28) 


9 22 






U3U3 


14 OO 




0.12326(47) 


9 19 


n (\7( ^fC\ 

U.O ( ^OD ) 


fl Q1 90^1^ 1 Af\~\ 
u.oizyo^i4D ) 


S\1L\ 


14 35 


1 55 ( 58^1 


0.24744(23) 


9 31 


fil f^fi 1 ! 

U . U 1 1 'J\J 1 


Sfi205f 65 s ! 


<j 1 u>2 


14 35 


1 54/581 


0.23160(25) 


9 28 


51 f35l 


35402 ( 741 


Si ti3 


14 35 


1.48(57) 


0.22326(29) 


9 27 


53(361 


35022( 81 1 


S2 til 


14 35 


1 53(581 


0.22251(24) 


9 29 


56(361 


34873( 721 




14 35 


1.50(57) 


0.20491(25) 


9 27 


0.41(32) 


34074( 81 1 


O J U> J, 


14 35 


1 43(561 


0.19546(28) 


9 25 


0.54(38) 


33699( 891 








(3 = 6.25 










Pseudo-Scalars 






Vectors 






range 




mass 


range 


2 / AT 

X /Ndf 


mass 


SlSl 


16 25 


1.62(91) 


0.24722(35) 


11 42 


1.24(43) 


0.32526( 48) 


S2S2 


16 25 


1.66(92) 


0.21638(36) 


11 39 


1.40(48) 


0.30620( 57) 


til til 


ID 20 


1. ( 4^yo j 


0.17228(37) 


11 33 


1. 04^00 ) 


fl oooon/ 1 A\ 


ti2 ti2 


1 £ op; 

ID ZO 


l.o4^y ( J 


0.13441(38) 


11 26 


1 Q r >f7f<i 


9fi7Q8i' 1 09^ 




1 £ OP» 
ID ZO 


1.0 / 1 y / ) 


0.10684(42) 


11 21 




9^09/1 f' 1 /I 


<5 1 U'l 


16 25 


1 fi8fQ^ 


0.21253(36) 


11 38 


1.33(47) 


S04'W 5Q^1 


■J 1 l*J 


16 25 


1.71(93) 


0.19806(38) 


11 35 


1.37(51) 


29683( 671 


Si ti3 


16 25 


1 72(931 


0.18943(40) 


11 32 


1 SI 


29281 ( 741 


S2 til 


16 25 


1 70(931 


0.19542(36) 


11 36 


1 48(521 


29467( 651 




16 25 


1.74(94) 


0.17976(38) 


11 33 


1.45(54) 


2871 1 ( 731 


S2 ti3 


16 25 


1 76(951 

. 1 VJ \ >J<J 1 


0.17032(40) 


11 30 


1 58(611 


28323( 81 1 








(3 = 6.47 










Pseudo-bcalars 






Vectors 






range 


2 / 7v r 

X /Ndf 


mass 


range 


2 / a r 


mass 


S1S1 


20 35 


0.37(38) 


0.18824(33) 


14 43 


1.37(58) 


0.24786( 55) 


S2S2 


20 35 


0.53(45) 


0.16314(35) 


14 37 


1.06(65) 


0.23029( 95) 


U1U1 


20 35 


0.88(57) 


0.12410(38) 


14 30 


1.07(58) 


0.21233( 91) 


U2U2 


20 35 


1.19(67) 


0.09988(41) 


14 26 


1.38(77) 


0.20265(119) 


M3M3 


20 35 


1.69(81) 


0.07678(45) 


14 21 


1.85(117) 


0.19656(186) 


SlWl 


20 35 


0.58(47) 


0.15892(37) 


14 35 


1.17(57) 


0.23012( 70) 


S1U2 


20 35 


0.69(50) 


0.14997(39) 


14 33 


1.02(54) 


0.22575( 78) 


SiU S 


20 35 


0.81(54) 


0.14307(43) 


14 30 


0.73(50) 


0.22336( 92) 


S2U1 


20 35 


0.68(51) 


0.14479(37) 


14 33 


1.20(58) 


0.22213( 76) 


S2U2 


20 35 


0.80(54) 


0.13502(39) 


14 31 


0.91(52) 


0.21794( 85) 


S2«3 


20 35 


0.93(57) 


0.12736(43) 


14 28 


0.91(59) 


0.21527(100) 



2G 



TABLE VII. The same as Table ^ for octet baryons. 



P = 5.90 







S-like 






A-like 






range 


X 2 /N df 


mass 


range 


X 2 /N df 


mass 


SlSlSl 


5 22 


0.47(35) 


0.7843(12) 








S2S2S2 


5 21 


0.48(35) 


0.7267(13) 








UlUlUl 


5 18 


0.41(36) 


0.6572(15) 








U2U2U2 


5 15 


0.44(44) 


0.6048(18) 








U3U3U3 


5 12 


0.28(46) 


0.5719(27) 








UlUlSl 


5 19 


0.35(33) 


0.7082(14) 


5 20 


0.51(38) 


0.6952(14) 


U2U2S1 


5 17 


0.34(35) 


0.6808(15) 


5 18 


0.55(42) 


0.6588(15) 


U3U3S1 


5 16 


0.26(31) 


0.6664(17) 


5 15 


0.70(48) 


0.6361(19) 


S1S1U1 


5 21 


0.49(36) 


0.7391(13) 


5 20 


0.40(34) 


0.7488(13) 


S1S1U2 


5 20 


0.63(41) 


0.7223(13) 


5 19 


0.37(34) 


0.7366(14) 


S1S1U3 


5 19 


0.76(47) 


0.7134(14) 


5 19 


0.45(37) 


0.7303(14) 


U1UIS2 


5 19 


0.40(35) 


0.6849(14) 


5 19 


0.51(39) 


0.6774(14) 


U2U2S2 


KIT 

5 17 


U. 00(00) 


U.dod4(1dJ 


KIT 


A E A/ A 0\ 

U.oU(4z) 


A C A AO/ 1 £i\ 

U.04Uz(l0J 


U3U3S2 


c 1 c 

lo 


U.ZO(OO) 


a a A 1 1 ( 1 o"\ 
U.0411(l5j 


K 1 A 
14 


A AO / A A\ 

U.4o(44) 


A CI TZ. 1 1 A\ 

u.d! /o(iy) 


S2S2U1 


5 20 


a a a 1 a \ 

0.45(36) 


A TA 1 "1 f "1 A \ 

0.7011(14) 


5 19 


0.37(34) 


A '7A'7C / 1 A \ 

0.7075(14) 


S2S2U2 


Kin 
19 


0.08(41 ) 


A COOI /'I A\ 

U.6831(14J 


K 1 A 
19 


A OO /O A \ 

0.30(34) 


A dCi A Q / A A\ 

0.0940(14) 


S2S2U3 


K 1 O 

lo 


A ro/)to\ 

0. 00(43 ) 


A t?TO d/lC^ 

0.6r34(loj 


K 1 O 

18 


0.40(30) 


A COOO/1 K \ 

0.OOO2(lO) 






E-like 


P = 6.10 




A-like 






range 


X 2 /N df 


mass 


range 


X 2 /N df 


mass 


S1S1S1 


7 28 


1.06(47) 


0.5892(10) 








S2S2S2 


7 26 


0.86(44) 


0.5383(11) 








UlUlUl 


7 23 


0.69(44) 


0.4831(12) 








U2U2U2 


7 19 


0.77(56) 


0.4421(15) 








U3U3U3 


7 16 


0.66(61) 


0.4179(19) 








U1U1S1 


7 24 


0.62(40) 


0.5255(11) 


7 25 


0.98(48) 


0.5150(11) 


U2U2S1 


7 21 


0.69(48) 


0.5043(13) 


7 23 


0.75(46) 


0.4855(13) 


U3U3S1 


7 19 


0.82(58) 


0.4933(14) 


7 20 


0.76(54) 


0.4685(15) 


S1S1U1 


7 26 


0.91(45) 


0.5515(10) 


7 26 


0.80(43) 


0.5594(11) 


S1S1U2 


7 25 


1.04(49) 


0.5384(11) 


7 25 


0.72(41) 


0.5497(11) 


S1S1U3 


7 25 


1.12(51) 


0.5315(11) 


7 24 


0.56(38) 


0.5448(12) 


M1W1S2 


7 24 


0.62(40) 


0.5050(12) 


7 24 


0.69(42) 


0.4989(12) 


U2U2S2 


T O 1 

r 21 


A P.f\l AQ\ 

0. 69(40 ) 


A /I OOO / 1 0\ 

U.4o2o(13J 


TOO 
( 22 


A TC / A o\ 

0. /D(40j 


a a an a / a o\ 
0.4094(13) 


U3U3S2 


Tin 
( 19 


A TO / K T*\ 

0. 10(01 ) 


A /I T1 O / 1 K\ 
0.4/ 12(10 ) 


T 1 A 
( 19 


A 71 /rr\ 

0. / 1(00) 


0.4o23(lO) 


S2S2U1 


t k 
( 20 


A AK / A T\ 
U.90(4/ ) 


A CI 01 / I 1 \ 

u. 0101(11 ) 


T O K 

( lo 


A QO/ A 9\ 

U.o2(4o ) 


A KOQ1 /1 1 \ 

U.02ol(ll ) 


S2S2U2 


T O /I 

f 24 


A TO i A A\ 

0. /3(44) 


A K AO /I / 1 0\ 

0.0034(12 ) 


T O A 

1 24 


A KO/OA\ 

0. 00(39 ) 


A K 1 OO / 1 0\ 

0.0120(12) 


S2S2U3 


TOO 

I 23 


A TO / A d\ 

0. /0(4t>) 


A /I AKT/ 1 0\ 

0.490/ (12J 


TOO 

( 23 


A K A /OA\ 

0. 04(39 ) 


A K ATA / 1 0^ 

0.00/9(12) 






E-like 


/3 = 6.25 




A-like 






range 


X 2 /N df 


mass 


range 


X 2 /N df 


mass 


S1S1S1 


8 34 


1.54(52) 


0.4920( 8) 








S2S2S2 


8 32 


1.45(55) 


0.4553( 9) 








UlUlUl 


8 28 


1.71(63) 


0.4069(11) 








U2U2U2 


8 24 


1.57(67) 


0.3694(13) 








U3U3U3 


8 20 


1.14(67) 


0.3451(16) 








UlUlSl 


8 29 


1.65(61) 


0.4407(10) 


8 30 


1.64(61) 


0.4319(10) 


U2U2S1 


8 26 


1.50(62) 


0.4208(12) 


8 28 


1.59(62) 


0.4050(11) 


U3U3S1 


8 24 


1.50(67) 


0.4093(13) 


8 25 


1.41(63) 


0.3875(12) 


S1S1U1 


8 32 


1.50(57) 


0.4611( 9) 


8 32 


1.41(53) 


0.4678( 9) 


S1S1U2 


8 31 


1.56(59) 


0.4485(10) 


8 30 


1.57(57) 


0.4588(10) 


S1S1U3 


8 30 


1.51(58) 


0.4408(10) 


8 28 


1.71(62) 


0.4538(10) 


U1U1S2 


8 29 


1.67(61) 


0.4262(10) 


8 29 


1.69(63) 


0.4207(10) 


u 2 u 2 s 2 


8 26 


1.50(62) 


0.4056(12) 


8 27 


1.63(64) 


0.3937(11) 


U3U3S2 


8 23 


1.57(71) 


0.3936(13) 


8 24 


1.44(66) 


0.3759(13) 


S2S2U1 


8 31 


1.55(58) 


0.4374(10) 


8 30 


1.58(59) 


0.4419(10) 
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S2S2W2 


zy 


l.UO^OO ) 




Q OS 
O ZO 


1. / o{Oo ) 


U.4ozo\ii ) 


S2S2W3 


ft OS 
Zo 


l.OD^Dz j 


n /i 1 ^7/ 1 1 ^ 

U.410 i l) 


ft 07 
Z 1 


1. / Z{DO ) 










13 = 6.47 












E-like 






A-like 






range 


9 (nr 


mass 


range 


9 1 HT 

X 2 /N df 


mass 


SlSlSl 


10 32 


0.92(60) 


0.3724(10) 








S2S2S2 


10 30 


0.93(62) 


0.3431(11) 








MlWlWl 


10 25 


1.19(82) 


0.3011(14) 








M2M2W2 


10 21 


1.18(86) 


0.2787(17) 








U3U3U3 


10 17 


1.33(101) 


0.2597(23) 








U1U1S1 


10 28 


0.97(66) 


0.3300(12) 


10 28 


1.07(71) 


0.3223(12) 


U2U2S1 


10 25 


1.06(71) 


0.3184(14) 


10 25 


1.25(84) 


0.3059(14) 


u 3 u :i si 


10 22 


1.04(72) 


0.3108(17) 


10 21 


1.49(94) 


0.2929(17) 


Si Si Ml 


10 30 


0.91(62) 


0.3471(11) 


10 30 


0.88(59) 


0.3525(11) 


S1S1M2 


10 29 


0.96(64) 


0.3397(11) 


10 29 


0.90(58) 


0.3473(12) 


S 1 S 1 U 3 


10 28 


0.99(65) 


0.3340(12) 


10 27 


0.93(59) 


0.3440(13) 


U1U1S2 


10 27 


1.05(71) 


0.3181(13) 


10 27 


1.12(76) 


0.3131(13) 


U2U2S2 


10 24 


1.14(76) 


0.3060(15) 


10 24 


1.31(87) 


0.2965(15) 


U3U3S2 


10 21 


1.00(77) 


0.2978(18) 


10 20 


1.45(97) 


0.2837(18) 


S2S2M1 


10 29 


1.00(67) 


0.3277(12) 


10 29 


0.96(64) 


0.3316(12) 


S2S2U2 


10 27 


1.10(75) 


0.3195(13) 


10 27 


1.02(66) 


0.3262(13) 


S2S2U3 


10 26 


1.11(79) 


0.3131(14) 


10 25 


1.04(67) 


0.3227(14) 
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TABLE VIII. The same as Table for decuplet baryons. 







P = 5.90 






= 6.10 






range 


X 2 /N df 


mass 


range 


X 2 /N df 


mass 


SlSlSl 


5 18 


0.41(37) 


0.8613(17) 


7 23 


0.83(49) 


0.6483(13) 


S2S2S2 


5 16 


0.38(39) 


0.8173(18) 


7 21 


0.71(48) 


0.6096(14) 


IlllilMl 


5 14 


0.59(55) 


0.7715(22) 


7 18 


1.12(69) 


0.5746(17) 


U2U2U2 


5 12 


0.96(81) 


0.7423(25) 


7 15 


1.12(80) 


0.5520(20) 


U3U3U3 


5 11 


0.99(92) 


0.7281(28) 


7 14 


1.11(86) 


0.5408(23) 


UlUlSl 


5 15 


0.40(43) 


0.8003(20) 


7 19 


0.80(56) 


0.5983(16) 


U2U2S1 


5 14 


0.48(50) 


0.7805(22) 


7 18 


1.13(70) 


0.5830(17) 




5 13 


0.48(53) 


0.7709(23) 


7 17 


1.25(77) 


0.5756(19) 


S1S1M1 


5 17 


0.35(36) 


0.8308(18) 


7 21 


0.73(49) 


0.6229(14) 


S1S1U2 


5 16 


0.35(37) 


0.8208(19) 


7 20 


0.74(51) 


0.6153(15) 


S1S1U3 


5 15 


0.34(39) 


0.8157(20) 


7 19 


0.78(55) 


0.6116(16) 


U1U1S2 


5 15 


0.49(48) 


0.7856(21) 


7 19 


0.87(58) 


0.5853(16) 


U2U2S2 


CIA 

5 14 


0.62(57) 


0.7658(23) 


7 17 


1.23(76) 


A f t~7f\ 1 /1 0\ 

0.5701(18) 


U3U3S2 


5 13 


0.61(60) 


0.7561(24) 


7 16 


A O O / C A \ 

0.83(64) 


0.5627(19) 


S2S2U1 


5 16 


a A 1 a o\ 

0.43(42) 


a orii A ( 1 C\\ 

0.8014(19) 


7 20 


0.75(52) 


0.5974(15) 


S2S2U2 


5 15 


A A A 1 A TL \ 

0.44(45) 


0.7912(20) 


7 19 


0.85(57) 


0.5897(16) 


S2S2U3 


C 1 A 

5 14 


A A Ct ( A C\\ 

0.46(49) 


0.7862(21) 


7 18 


A A^7 1 Ct A \ 

0.97(64) 


0.5859(17) 






P = 6.25 






P = 6.47 






range 


X 2 /N dS 


mass 


range 


X 2 /N df 


mass 


SiS 1 S 1 


8 27 


1.03(50) 


0.5405(11) 


10 29 


0.87(54) 


0.4133(15) 


S2S2S2 


8 25 


1.20(58) 


0.5124(12) 


10 27 


0.88(56) 


0.3913(17) 


U1U1U1 


8 21 


1.51(76) 


0.4795(15) 


10 22 


0.66(57) 


0.3639(21) 


U2U2U2 


8 19 


1.44(81) 


0.4573(18) 


10 20 


0.72(67) 


0.3514(24) 


U3U3U3 


8 17 


0.86(69) 


0.4455(20) 


10 17 


1.15(101) 


0.3435(29) 


UlUlSl 


8 23 


1.35(66) 


0.4991(13) 


10 26 


0.84(56) 


0.3804(18) 


U2U2S1 


8 21 


1.47(76) 


0.4842(15) 


10 24 


0.70(52) 


0.3715(20) 


U3U3S1 


8 20 


1.29(75) 


0.4760(16) 


10 22 


0.71(57) 


0.3667(22) 


S1S1M1 


8 25 


1.19(58) 


0.5198(12) 


10 28 


0.97(57) 


0.3970(17) 


S1S1W2 


8 24 


1.29(62) 


0.5123(13) 


10 27 


0.92(55) 


0.3932(17) 




8 23 


1 29(65") 


5080(13") 


10 26 


92(57) 


3906(18") 


U1U1S2 


8 22 


1.45(71) 


0.4897(14) 


10 25 


0.95(59) 


0.3730(19) 


U2U2S2 


8 21 


1.46(76) 


0.4747(15) 


10 23 


0.77(57) 


0.3640(21) 


U3U3S2 


8 19 


1.36(80) 


0.4666(17) 


10 21 


0.74(61) 


0.3590(23) 


S2S2U1 


8 24 


1.33(63) 


0.5010(13) 


10 26 


0.81(56) 


0.3818(18) 


S2S2U2 


8 23 


1.35(67) 


0.4933(14) 


10 25 


0.94(59) 


0.3780(18) 


S2S2U3 


8 22 


1.37(71) 


0.4892(14) 


10 24 


0.69(52) 


0.3746(20) 
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TABLE IX. Bare AWI quark masses 2m q obtained at the simulation points. Here, 2m q is a short cut notation for 

m q/ +m qg . The fit range [t min ,t max ] is also given. 







P = 5.90 




= 6.10 


P = 6.25 




P = 6.47 


K\K2 


range 


2mf wm 


range 


2m^ wl{0) 


range 2m q WI( - ^ 


range 


2mf wm 


s 1 s 1 


8 25 


0.090662( 82) 


10 32 


0.069701( 56) 


10 39 0.060678(48) 


10 53 


0.046245( 46) 


S2S2 


8 25 


0.067556( 77) 


10 32 


0.050123( 53) 


10 39 0.046676(47) 


10 53 


0.034949( 44) 


MlWl 


8 25 


0.042373( 72) 


10 32 


0.030977( 50) 


10 39 0.029459(44) 


10 53 


0.020087( 41) 


U2U2 


8 25 


0.025840( 74) 


10 32 


0.018363( 50) 


10 39 0.017590(42) 


10 53 


0.012752( 39) 


U3U3 


8 25 


0.016953(201) 


10 32 


0.011828( 68) 


10 39 0.010777(44) 


10 53 


0.007253( 40) 


SlUl 


8 25 


0.066309( 78) 


10 32 


0.050203( 54) 


10 39 0.044982(47) 


10 53 


0.033101( 45) 


S1W2 


8 25 


0.057951( 78) 


10 32 


0.043825( 54) 


10 39 0.038975(46) 


10 53 


0.029392( 44) 


S1U3 


8 25 


0.053722( 82) 


10 32 


0.040611( 56) 


10 39 0.035554(47) 


10 53 


0.026619( 45) 


S2U1 


8 25 


0.054911( 75) 


10 32 


0.040518( 52) 


10 39 0.038040(46) 


10 53 


0.027495( 43) 




2o 


U.04bt)ll( to) 


1 A OO 

10 61 


0.034180( 52) 


10 39 0.032063(45) 


10 53 


0.023802( 43J 


S2U3 


8 25 


0.042394( 79) 


10 32 


0.030978( 54) 


10 39 0.028654(45) 


10 53 


0.021038( 43) 






TABLE X 


. Values of x, y 


and yf used in the ratio tests for m 2 PS and fps- 












13 = 5.90 






K1K2 






X 




V 




Vf 


SlWl 




-0.09551( 24) 




0.98835( 38) 




1.00408( 59) 


S1U2 




-0.25594( 87) 




0.9749 (12 ) 




1.0122 (17 ) 


S1M3 




-0.4480 (59 ) 




0.9558 (62 ) 




1.0313 (80 ) 






-0.03613( 10) 




0.99602( 17) 




1.00166( 24) 


S2M2 




-0.15161( 60) 




0.98710( 84) 




1.0076 (12 ) 


S2W3 




-0.3089 (50 ) 




0.9725 (57 ) 




1.0249 (75 ) 










P = 6.10 






K1K2 






X 








Vf 


SlWl 




-0.10843( 25) 




0.98764( 58) 




1.00369( 82) 


S1U2 




-0.28809( 90) 




0.9718 (18 ) 




1.0128 (25 ) 


S1W3 




-0.4988 (29 ) 




0.9553 (62 ) 




1.0346 (75 ) 


S2W1 




-0.03845( 10) 




0.99589( 24) 




1.00158( 31) 


S2U2 




-0.16528( 61) 




0.9852 (13 ) 




1.0086 (17 ) 


S2W3 




-0.3360 (23 ) 




0.9725 (58 ) 




1.0284 (66 ) 










P = 6.25 






Kl«2 






X 








Vf 


SlWl 




-0.08628( 19) 




0.98986( 53) 




1.00232( 70) 


S1U2 




-0.24925( 69) 




0.9712 (16 ) 




1.0094 (21 ) 


S1M3 




-0.4746 (19 ) 




0.9456 (34 ) 




1.0243 (47 ) 






-0.03518( 10) 




0.99583( 24) 




1.00116( 30) 


S 2 U 2 




-0.15627( 49) 




0.9817 (12 ) 




1.0068 (14 ) 


S2U3 




-0.3459 (15 ) 




0.9599 (29 ) 




1.0201 (38 ) 










P = 6.47 






n 1 K 2 






X 








Vf 


SlWl 




-0.11458( 32) 




0.9870 (15 ) 




1.0054 (17 ) 


S1U2 




-0.26924( 95) 




0.9699 (37 ) 




1.0145 (43 ) 


Sl«3 




-0.5417 (27 ) 




0.9402 (81 ) 




1.0336 (99 ) 


S2U1 




-0.05086( 16) 




0.99414( 75) 




1.00272( 83) 


S2U2 




-0.16660( 67) 




0.9808 (26 ) 




1.0096 (30 ) 


S2U3 




-0.3961 (23 ) 




0.9546 (69 ) 




1.0262 (82 ) 
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TABLE XI. Values of S and a* from the fit (|l|). 



p 


One parameter fit 
y — 1 + Sx 

S X 2 /N df 


S 


Two parameter fit 
y = 1 + Sx + axz 
ctx Q$ 


X 2 /N df 


5.90 


0.106( 5) 


7.0 


0.016(12) 


2.69(29) 


0.59( 6) 


0.61 


6.10 


0.103( 6) 


2.2 


0.042(21) 


1.88(55) 


0.40(12) 


0.12 


6.25 


0.117( 7) 


0.1 


0.112(15) 


0.14(36) 


0.03( 8) 


0.03 


6.47 


0.113(13) 


0.0 


0.113(33) 


0.01(88) 


0.00(19) 


0.02 



TABLE XII. Parameters from Q^PT chiral fit (|28[) with q$ = for pseudoscalar meson masses. 









degenerate fit 








P = 5.90 




P = 6.10 


P = 6.25 


P = 6.47 


ftc 


0.1598315(68) 




0.1549903(55) 


0.1525355(69) 


0.1498049(79) 


A [lattice] 


1.1127(39) 




0.8784(42) 


0.7396(67) 


0.5862(85) 


A [GeV] 


2.152(19) 




2.231(22) 


2.271(31) 


2.322(59) 


6 


0.1061(45) 




0.0815(68) 


0.0598(96) 


0.079(20) 


B 


1.012(29) 




1.096(47) 


1.088(55) 


1.42(15) 


A x [GeV] 


0.761(7) 




0.789(8) 


0.803(11) 


0.821(21) 






non-de. 


generate fit with n\ — si 








P = 5.90 




P = 6.10 


P = 6.25 


P = 6.47 


A [lattice] 


1.1161(49) 




0.8740(49) 


0.7339(62) 


0.5771(90) 


A [GeV] 


2.158(20) 




2.220(23) 


2.254(30) 


2.286(58) 


S 


0.0943(57) 




0.0813(62) 


0.0774(88) 


0.090(24) 


B 


0.893(51) 




1.121(55) 


1.285(66) 


1.62(21) 






non-dej 


generate fit with «i = S2 








P = 5.90 




P = 6.10 


P = 6.25 


P = 6.47 


A [lattice] 


1.1177(54) 




0.8731(55) 


0.7312(65) 


0.5731(96) 


A [GeV] 


2.161(20) 




2.218(23) 


2.246(30) 


2.270(59) 


S 


0.0980(63) 




0.0827(64) 


0.0815(99) 


0.101(25) 


B 


0.906(66) 




1.143(65) 


1.360(79) 


1.78(24) 


TABLE XIII. 


Values of S obtained from ratio tests for mps and fps and Q^PT chiral fits. DG denotes degenerate fit and 


ND,Si non-degenerate fits with one of k 


being fixed to 


Si. 










P = 5.90 


P = 6.10 


P = 6.25 


P = 6.47 


ratio 












m PS 




0.106( 5) 


0.103( 6) 


0.117( 7) 


0.113(13) 


fps 




0.095(13) 


0.091(17) 


0.079(17) 


0.111(32) 


fit 


2^7 A£ 










DG 


4 


0.122( 4) 


0.091( 5) 


0.065(11) 


0.088( 9) 


DG 


8 


0.114( 5) 


0.086( 8) 


0.062(10) 


0.083(23) 


DG 


16 


0.106( 5) 


0.082( 7) 


0.060(10) 


0.079(20) 


DG 


32 


0.100( 4) 


0.077( 6) 


0.058( 9) 


0.075(19) 


ND,si 


16 


0.094( 6) 


0.081( 6) 


0.077( 9) 


0.090(24) 


ND,s 2 


16 


0.098( 6) 


0.083( 6) 


0.082(10) 


0.101(25) 
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TABLE XIV. Parameters from the QxPT fit (^) for vector meson masses. 



X 2 /N df 
m v (Lat) 
m° v (GeV) 
Ci (Lat) 
Ci (GeV" 1 ) 

Cl/2 


ft k on 

p=o.aU 
0.40 
0.3969(49) 
0.7675(32) 
0.919(38) 
0.475(23) 
-0.058(27) 


ft a i n 
p=6.10 

0.07 

0.3030(38) 

0.7695(33) 

1.259(54) 

0.495(25) 

-0.075(28) 


ft a oc 
p=6.2o 

0.08 
0.2502(40) 
0.7684(42) 
1.496(79) 
0.487(30) 
-0.065(35) 


ft HAT 

p=0.4( 
0.72 
0.1966(59) 
0.7789(85) 
2.27(22) 
0.574(66) 
-0.155(72) 


TABLE XV. Parameters from the QxPT fit, @ and (|l|) for octet baryon 
coefficient of the mps-linear term in the degenerate mass formula. 


masses. C1/2 = ~(3nS/2)(D — 


3F) 2 is the 


X 2 /df 
m (Lat) 
m (GeV) 
bp (Lat) 
bp (GeV -1 ) 
bo (Lat) 
b D (GeV" 1 ) 
F 
D 

Cl/2 


/3=5.90 
0.42 
0.5141(29) 
0.9941(84) 
-0.534(27) 
-0.276(14) 
0.034(13) 
0.018(7) 
0.334(14) 
0.484(13) 
-0.126(15) 


(3=6.10 
0.60 
0.3774(25) 
0.9585(81) 
-0.696(36) 
-0.274(14) 

U.Uo2(i / ) 

0.020(7) 
0.326(14) 
0.475(15) 
-0.120(15) 


(3=6.25 
0.73 
0.3075(20) 
0.945(10) 
— U.o/o(o / ) 
-0.268(13) 
0.071(20) 
0.023(6) 
0.315(13) 
0.450(13) 
-0.116(13) 


(3=6.47 
0.10 

0.2353(33) 
0.932(19) 
— u.y 1 

-0.245(30) 
0.159(60) 
0.040(15) 
0.299(31) 
0.467(28) 

-0.087(30) 


TABLE XVI. Parameters from the QxPT fit (fi 
mps-linear term in the degenerate mass formula. 


^) for decuplet baryon masses. 


C1/2 = — (5tt8/6)H 2 is the coefficient of the 


X 2 /df 
m° D (Lat) 
m° D (GeV) 
c (Lat) 
c (GOV" 1 ) 
H 2 
C 2 

Ci/2 


/3=5.90 
0.03 
0.7144(64) 
1.381(15) 
0.483(16) 
0.250(9) 
0.65(13) 
-0.49(26) 
-0.169(33) 


/3=6.10 
0.03 
0.5301(52) 
1.346(14) 
0.630(27) 
0.248(11) 
0.55(15) 
-0.42(30) 
-0.145(39) 


/3=6.25 
0.03 
0.4309(48) 
1.323(17) 
0.735(31) 
0.239(11) 
0.39(15) 
-0.75(28) 
-0.101(40) 


(3=6.47 
0.04 
0.3390(87) 
1.343(36) 
1.028(110) 
0.259(29) 
0.70(40) 
0.88(79) 
-0.18(10) 
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TABLE XVII. Hadron masses in lattice units from QxPT fits. 







K input 








(3=5.90 


/3=6.10 


/3=6.25 


(3=6.47 




0.1598315( 68) 


0.1549903( 55) 


0.1525355( 69) 


0.1498049( 79) 


Kud 


0.1597449( 63) 


0.1549274( 47) 


0.1524834( 57) 


0.1497687( 65) 


K s 


0.157104 (49 ) 


0.153103 (33 ) 


0.151056 (33 ) 


0.148730 (49 ) 


7T 


0.06982( 58) 


0.05315( 45) 


0.04396( 48) 


0.03408( 68) 


P 


0.3974 (33 ) 


0.3025 (26 ) 


0.2502 (27 ) 


0.1940 (38 ) 


K 


0.2574 (22 ) 


0.1959 (17 ) 


0.1621 (18 ) 


0.1256 (25 ) 


K* 


0.4427 (23 ) 


0.3376 (18 ) 


0.2797 (19 ) 


0.2150 (27 ) 


Vs 


0.3510 (32 ) 


0.2701 (24 ) 


0.2240 (26 ) 


0.1740 (41 ) 


<f> 


0.4899 (21 ) 


0.3745 (17 ) 


0.3106 (18 ) 


0.2384 (27 ) 


N 


0.5202 (32 ) 


0.3824 (27 ) 


0.3117 (21 ) 


0.2391 (35 ) 


A 


0.7096 (46 ) 


0.5277 (37 ) 


0.4308 (35 ) 


0.3363 (58 ) 


A 


0.5840 (30 ) 


0.4340 (24 ) 


0.3543 (21 ) 


0.2713 (33 ) 


S 


0.6294 (31 ) 


0.4691 (27 ) 


0.3841 (25 ) 


0.2948 (38 ) 


s* 


0.7469 (36 ) 


0.5588 (30 ) 


0.4576 (27 ) 


0.3578 (41 ) 


A ss 


0.6986 (35 ) 


0.5241 (30 ) 


0.4308 (30 ) 


0.3295 (45 ) 




0.6671 (35 ) 


0.5001 (28 ) 


0.4091 (27 ) 


0.3140 (44 ) 


H* 


0.7881 (34 ) 


0.5924 (29 ) 


0.4870 (26 ) 


0.3791 (40 ) 


N sss 


0.7486 (40 ) 


0.5641 (32 ) 


0.4640 (34 ) 


0.3554 (52 ) 


Q 


0.8333 (36 ) 


0.6287 (29 ) 


0.5189 (29 ) 


0.4003 (45 ) 






4> input 








/3=5.90 


13=6.10 


/3=6.25 


(3=6.47 


Kg 


0.156156 (105 ) 


0.152519 ( 69 ) 


0.150621 ( 70 ) 


0.148393 ( 86 ) 


K 


0.2965 (41 ) 


0.2232 (31 ) 


0.1836 (34 ) 


0.1433 (43 ) 


K* 


0.4611 (34 ) 


0.3508 (26 ) 


0.2902 (28 ) 


0.2241 (38 ) 


Vs 


0.4092 (61 ) 


0.3110 (46 ) 


0.2567 (50 ) 


0.2011 (66 ) 




0.5272 (44 ) 


0.4013 (34 ) 


0.3319 (36 ) 


0.2573 (51 ) 


A 


0.6094 (40 ) 


0.4525 (31 ) 


0.3689 (30 ) 


0.2833 (43 ) 


E 


0.6620 (45 ) 


0.4930 (38 ) 


0.4038 (39 ) 


0.3109 (53 ) 


S* 


0.7639 (42 ) 


0.5711 (35 ) 


0.4675 (32 ) 


0.3666 (48 ) 


A ss 


0.7490 (60 ) 


0.5612 (47 ) 


0.4618 (53 ) 


0.3542 (67 ) 




0.7178 (61 ) 


0.5372 (48 ) 


0.4387 (51 ) 


0.3391 (71 ) 




0.8232 (52 ) 


0.6177 (41 ) 


0.5073 (42 ) 


0.3968 (58 ) 


N sss 


0.8134 (72 ) 


0.6118 (56 ) 


0.5034 (63 ) 


0.3878 (83 ) 


Q 


0.8876 (67 ) 


0.6677 (51 ) 


0.5502 (55 ) 


0.4267 (75 ) 
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TABLE XVIII. Hadron masses in units of GeV from QxPT fits. 



K input 





exp. 


/3=5.90 


/3=6.10 


/3=6.25 


/3=6.47 


a' 1 




1.934(16) 


2.540(22) 


3.071(34) 


3.961(79) 


K* 


0.8961 


0.8561 ( 29) 


0.8575( 30) 


0.8591 ( 40) 


0.8516( 71) 






0.6786( 13) 


0.6860( 14) 


0.6879( 19) 


0.6893( 59) 


<t> 


1.0194 


0.9472( 43) 


0.9512( 44) 


0.9538( 58) 


0.9445(102) 


N 


0.9396 


1.0060( 87) 


0.9712( 83) 


0.9574(106) 


0.9472(194) 


A 


1.2320 


1.3722(128) 


1.3404(116) 


1.3229(142) 


1.3321(278) 


A 


1.1157 


1.1294( 73) 


1.1023( 69) 


1.0881( 88) 


1.0746(150) 


£ 


1.1926 


1.2171( 65) 


1.1915( 66) 


1.1797( 82) 


1.1676(138) 


S* 


1.3837 


1.4442(104) 


1.4193( 97) 


1.4055(120) 


1.4173(215) 


A ss 




1.3509( 60) 


1.3312( 64) 


1.3232( 75) 


1.3052(132) 




1.3149 


1.2899( 60) 


1.2702( 59) 


1.2563( 75) 


1.2440(122) 


5* 


1.5318 


1.5240( 90) 


1.5048( 89) 


1.4957(111) 


1.5019(194) 


N sss 




1.4475( 56) 


1.4327( 58) 


1.4250( 70) 


1.4077(120) 


n 


1.6725 


1.6113( 82) 


1.5970( 84) 


1.5935(106) 


1.5858(185) 








4> input 








exp. 


/3=5.90 


/3=6.10 


/3=6.25 


/3=6.47 


K 


0.4977 


0.5733( 34) 


0.5669( 34) 


0.5640( 46) 


0.5675( 72) 


K* 


0.8961 


0.8917( 11) 


0.8910( 11) 


0.8914( 14) 


0.8878( 29) 


r)s 




0.7913( 54) 


0.7898( 52) 


0.7883( 71) 


0.7967(109) 


N 


0.9396 


1.0060( 87) 


0.9712( 83) 


0.9574(106) 


0.9472(194) 


A 


1.2320 


1.3722(128) 


1.3404(116) 


1.3229(142) 


1.3321(278) 


A 


1.1157 


1.1784( 48) 


1.1492( 46) 


1.1331( 55) 


1.1221(102) 


£ 


1.1926 


1.2800( 39) 


1.2521( 44) 


1.2402( 47) 


1.2317( 85) 


£* 


1.3837 


1.4771( 87) 


1.4506( 81) 


1.4357( 98) 


1.4523(183) 


A ss 




1.4484( 28) 


1.4254( 34) 


1.4183( 33) 


1.4032( 65) 




1.3149 


1.3879( 26) 


1.3644( 28) 


1.3473( 31) 


1.3432( 55) 




1.5318 


1.5918( 56) 


1.5690( 57) 


1.5580( 65) 


1.5717(128) 


N sss 




1.5729( 27) 


1.5540( 26) 


1.5461( 38) 


1.5362( 61) 




1.6725 


1.7162( 33) 


1.6959( 33) 


1.6897( 38) 


1.6903( 93) 



TABLE XIX. Coefficients S in units of GeV of O(a) terms in linear continuum extrapolations of hadron masses, x /Ndf of 
the fits are also shown. 



input 




K 


K* 





N 


A 


S 




A 


E* 




n 


m K 


S 




-0.004 


-0.018 


0.278 


0.208 


0.173 


0.144 


0.173 


0.118 


0.083 


0.062 




X 2 /N df 




0.47 


0.45 


0.15 


0.01 


0.03 


0.03 


0.36 


0.52 


0.28 


0.03 




S 


0.067 


0.007 




0.278 


0.214 


0.169 


0.149 


0.173 


0.120 


0.092 


0.079 






0.35 


0.49 




0.15 


0.13 


0.38 


1.29 


0.36 


1.08 


1.63 


1.02 
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TABLE XX. Comparison of statistical errors (Stat.), estimated 0(a 2 ) effects of continuum extrapolations and estimated 
systematic errors of chiral extrapolations (mass difference from QxPT and polynomial chiral fits) . Deviations of the experimental 
value from our final mass prediction are also given. Errors or deviations quoted in percent are relative to the central values of 
the predicted mass; those in terms of a are normalized by statistical errors. 



K input 





Stat. 


0{a 2 ) error 


chiral fit error 


deviation 


K* 


1.04% 


0.00% 0.00a 


0.67% 0.65cr 


4.42% 4.26cr 


4> 


1.35% 


0.00% 0.00a 


1.35% l.OOcr 


6.54% 4.83cr 


N 


2.84% 


1.12% 0.39a 


1.05% 0.37a 


7.01% 2.47a 


A 


2.00% 


0.63% 0.31a 


0.24% 0.12cr 


9.46% 4.74a 


E 


1.67% 


0.43% 0.26a 


-0.63% 0.38cr 


6.80% 4.07a 




1.40% 


0.30% 0.21a 


0.88% 0.63a 


9.49% 6.78a 


A 


2.80% 


0.43% 0.15a 


-0.38% 0.14cr 


-2.00% 0.72a 


E* 


2.10% 


0.20% 0.10a 


-0.54% 0.26a 


1.84% 0.88a 


■H * 


1.76% 


0.10% 0.06(7 


0.29% 0.16cr 


4.97% 2.83a 


Q 


1.54% 


0.06% 0.04ct 


-0.07% 0.04a 


7.17% 4.65a 


4> input 




Stat. 


O{o i ) error 


chiral fit error 


deviation 


K 


1.78% 


0.06% 0.04(7 


-1.34% 0.75cr 


-10.06% 5.64a 


K* 


0.37% 


0.00% 0.00(7 


0.23% 0.63a 


0.84% 2.27a 


A 


1.27% 


0.66% 0.52a 


-0.21% 0.17cr 


5.21% 4.11a 


E 


0.96% 


0.41% 0.43a 


-1.54% 1.61cr 


1.41% 1.48a 




0.58% 


0.32% 0.55a 


-0.23% 0.40cr 


2.06% 3.52a 


E* 


1.72% 


0.21% 0.12a 


-0.53% 0.31cr 


-0.28% 0.17a 




1.05% 


0.12% 0.12a 


0.09% 0.08cr 


0.99% 0.94a 


Q 


0.60% 


0.09% 0.15a 


-0.64% 1.07cr 


1.53% 2.56a 
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TABLE XXI. Hadron masses in units of GeV from polynomial fits. 



K input 





exp. 


/3=5.90 


/3=6.10 


/3=6.25 


/3=6.47 


a' 1 




1.970( 8) 


2.584(13) 


3.109(20) 


4.102(60) 


K* 


0.8961 


0.8726( 20) 


0.8722( 20) 


0.8678( 33) 


0.8672( 73) 






0.6906( 12) 


0.6930( 12) 


0.6915( 17) 


0.6908( 59) 


<t> 


1.0194 


0.9659( 19) 


0.9672( 22) 


0.9676( 25) 


0.9657( 65) 


N 


0.9396 


1.0139( 64) 


0.9833( 73) 


0.9640( 97) 


0.9576(196) 


A 


1.2320 


1.3968( 76) 


1.3621( 83) 


1.3365( 92) 


1.3516(241) 


A 


1.1157 


1.1600( 47) 


1.1265( 55) 


1.1045( 71) 


1.1002(139) 


£ 


1.1926 


1.2272( 44) 


1.1982( 50) 


1.1804( 61) 


1.1787(135) 


S* 


1.3837 


1.4626( 64) 


1.4342( 68) 


1.4160( 77) 


1.4311(201) 


A ss 




1.3629( 41) 


1.3384( 45) 


1.3245( 52) 


1.3249(127) 




1.3149 


1.3203( 42) 


1.2971( 45) 


1.2759( 53) 


1.2752(136) 


5* 


1.5318 


1.5542( 56) 


1.5290( 62) 


1.5169( 73) 


1.5310(190) 


N sss 




1.4790( 37) 


1.4608( 43) 


1.4462( 49) 


1.4348(133) 


n 


1.6725 


1.6451( 55) 


1.6262( 58) 


1.6098( 67) 


1.6118(166) 








4> input 








exp. 


/3=5.90 


/3=6.10 


/3=6.25 


/3=6.47 


K 


0.4977 


0.5572( 20) 


0.5541( 22) 


0.5527( 26) 


0.5534( 65) 


K* 


0.8961 


0.8998( 13) 


0.8989( 14) 


0.8943( 27) 


0.8953( 59) 


T)s 




0.7787( 28) 


0.7775( 32) 


0.7750( 35) 


0.7762( 86) 


N 


0.9396 


1.0139( 64) 


0.9833( 73) 


0.9640( 97) 


0.9576(196) 


A 


1.2320 


1.3968( 76) 


1.3621( 83) 


1.3365( 92) 


1.3516(241) 


A 


1.1157 


1.1969( 40) 


1.1638( 48) 


1.1413( 62) 


1.1383(132) 


£ 


1.1926 


1.2738( 35) 


1.2455( 40) 


1.2275( 54) 


1.2245(110) 


£* 


1.3837 


1.4881( 55) 


1.4615( 59) 


1.4426( 68) 


1.4566(175) 


A ss 




1.4367( 29) 


1.4131( 33) 


1.3982( 44) 


1.4011(104) 




1.3149 


1.3926( 31) 


1.3707( 32) 


1.3478( 45) 


1.3483(124) 




1.5318 


1.6035( 43) 


1.5814( 47) 


1.5677( 55) 


1.5848(156) 


N sss 




1.5805( 26) 


1.5591( 28) 


1.5453( 39) 


1.5355( 95) 


a 


1.6725 


1.7251( 36) 


1.7057( 35) 


1.6882( 45) 


1.6902(107) 



TABLE XXII. Hadron spectrum from linear continuum extrapolations of masses determined by polynomial chiral fits. 





rriK input 

mass(GeV) 


X 2 /N df 


input 

mass (GeV) 


X 2 /N df 


K 






0.546(06) 


0.07 


K* 


0.864(07) 


0.42 


0.891(05) 


0.70 




0.970(06) 


0.08 






N 


0.887(22) 


0.18 


0.887(22) 


0.18 


A 


1.022(16) 


0.60 


1.058(14) 


0.80 


S 


1.110(14) 


0.62 


1.158(12) 


0.74 




1.212(13) 


0.89 


1.285(11) 


1.73 


A 


1.252(24) 


1.07 


1.252(24) 


1.07 


S* 


1.351(20) 


1.14 


1.380(17) 


1.34 




1.463(18) 


1.05 


1.518(14) 


1.71 


Q 


1.560(17) 


0.43 


1.637(11) 


1.14 
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TABLE XXIII. Values of k awi from the quark mass based on the axial- vector Ward id entity and k c from various chiral 
extrapolations of pseudoscalar meson masses. The second column shows fit number in Table XLVI. 



(3 = 5.90 = 6.10 /3 = 6.25 = 6.47 

k awi (linear) 0.1598254(38) 0.1549765(21) 0.1525172(13) 0.1497968(11) 

K AWI (quadratic) 0.1598490(93) 0.1549821(23) 0.1525226(15) 0.1498001(12) 

^c(QxPT) Fit 5 0.1598315(68) 0.1549903(55) 0.1525355(69) 0.1498049(79) 

Kc(quadratic) Fit 2 0.1599531(73) 0.1550430(50) 0.1525650(53) 0.1498244(78) 

K c (cubic) Fit 3 0.1598993(63) 0.1550251(45) 0.1525579(67) 0.1498210(76) 



TABLE XXIV. Parameters of chiral fits to quark masses based on the axial-vector Ward identity. 





0=5.90 


0=6.10 


0=6.25 


0=6.47 


A dg 

A dg 
K AWI 

r>dg 
°2 


-4.3885(39) 
0.70142(62) 
0.1598299(31) 
0.6986(49) 
0.022(28) 


-4.8871(38) 
0.75739(59) 
0.1549768(18) 
0.7386(70) 
0.163(58) 


-5.1724(35) 
0.78888(53) 
0.1525160(13) 
0.753(12) 
0.37(12) 


-5.5069(47) 
0.82492(71) 
0.1497959(11) 
0.799(21) 
0.39(30) 


A S1 
A o 

A? 


-2.1705(16) 
0.35410(26) 


-2.4283(17) 
0.38170(26) 


-2.5771(15) 
0.39765(23) 


-2.7510(21) 
0.41553(31) 


A 32 

A\ 2 


-2.1554(22) 
0.34989(35) 


-2.4146(22) 
0.37809(34) 


-2.5654(19) 
0.39481(28) 


-2.7399(27) 
0.41304(40) 



TABLE XXV. Light quark masses in lattice units. 



m VWI 
m ud 

m AWI 
m ud 

mY WI (K) 
mf WI (K) 



0=5.90 
0.001696(37) 
0.001185(23) 
0.0543(10) 
0.0379( 7) 
0.0736(21) 
0.0512(15) 



0=6.10 
0.001311(35) 
0.000969(21) 
0.0398( 7) 
0.0296( 5) 
0.0523(15) 
0.0389(11) 



0=6.25 
0.001120(42) 
0.000845(25) 
0.0321( 8) 
0.0246( 6) 
0.0417(15) 
0.0321(12) 



0=6.47 
0.000808(64) 
0.000646(42) 
0.0241(12) 
0.0195( 9) 
0.0318(20) 
0.0258(16) 



TABLE XXVI. Light quark masses in units of GeV in the MS scheme at fi = 2 GeV. 



m VWI 

m AWI 
m ud 

mY WI (K) 
mf WI (K) 
mY WI (<P) 



0=5.90 
0.00418( 6) 
0.00302( 4) 
0.1338(13) 
0.0965( 9) 
0.1814(35) 
0.1304(25) 



0=6.10 
0.00431 ( 9) 
0.00327( 5) 
0.1306(13) 
0.1000( 8) 
0.1716(32) 
0.1315(25) 



0=6.25 
0.00450(15) 
0.00348( 8) 
0.1289(17) 
0.1014(11) 
0.1673(41) 
0.1322(34) 



0=6.47 
0.00424(28) 
0.00347(18) 
0.1267(37) 
0.1049(26) 
0.1669(67) 
0.1385(54) 
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TABLE XXVII. 


Parameters of continuum extrapolations of quark masses. Parameters are defined by mY WI 


= m + A VWI a 


and mf WI = m + A Awl a. 










m (MeV) 


j^VWI 




X 2 /N df 


m ud 


4.29(14) 


-0.12(31) 


-2.49(30) 


1.13 


m s (K) 


113.8(2.3) 


41.2(5.5) 


34.6(5.3) 


0.81 


m s (cj>) 


142.3(5.8) 


76(15) 


-25(14) 


0.31 


TABLE XXVIII. Plaquette values 


/3 = 5.90 




,3 = 6.10 = 6.25 




(3 = 6.47 


0.58184394(109) 


0.60413086( 71) 0.61810507( 55) 




0.63608939( 47) 


TABLE XXIX. A^ in units of MeV. 


method 


(3 = 5.90 


/3 = 6.10 (3 = 6.25 


(3 = 6.47 


(3 = oo 


O.MF 


233.1(1.9) 


231.4(2.0) 229.4(2.5) 


223.2(4.4) 


220.4(5.7) 


ay 


219.1(1.8) 


220.4(1.9) 220.2(2.4) 


216.1(4.3) 


219.5(5.4) 


a P w/o BLM 


258.3(2.2) 


249.8(2.1) 244.9(2.7) 


235.5(4.7) 


219.3(6.2) 


a p w/ BLM 


254.6(2.1) 


246.7(2.1) 242.1(2.6) 


233.1(4.6) 


218.0(6.1) 


a p — ► Ap 


258.2(2.2) 


249.8(2.1) 244.9(2.7) 


235.5(4.7) 


219.2(6.2) 
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TABLE XXX. Pseudoscalar meson decay constant, Z K f PS , where Z K is the K-dependent part of the renormalization constant 
defined in Eq.(F3). Ranges for fits in Eq.(Fol) are also given. 





range 


P = 5.90 


value 




P = 6.10 

range 


value 


SlSj 

S2S2 

U2U2 
U3U3 
S1U1 
S1M2 
S1M3 
S 2 Ml 
S2«2 
S2W3 


12 28 
12 28 
12 28 
12 28 
12 28 
12 20 
12 28 
12 28 
12 28 
12 28 
12 28 




0.10128(51) 
0.09379(53) 
0.08498(55) 
0.07875(59) 
0.07469(87) 

0.0929/(53) 

0.08985(55) 
0.08833(56) 
0.08935(54) 
0.08627(55) 
0.08473(57) 




12 35 
12 35 
12 35 
12 35 
12 35 

12 35 

12 35 
12 35 
12 35 
12 35 
12 35 


0.07412(44) 
0.06834(47) 
0.06232(53) 
0.05816(61) 
0.05553(76) 

0.06809(45) 
0.06608(51) 

0.06526(55) 
0.06532(50) 
0.06332(53) 
0.06247(57) 




range 


P = 6.25 


value 




P = 6.47 

range 


value 


Si Si 

S2S2 

MlUl 
U2U2 
U3U3 
S1U1 

S1U2 
S1U3 
S2W1 

S2^2 
S2M3 


17 42 
17 42 
17 42 

17/19 

17 42 
17 42 
17 42 
17 42 
17 42 
17 42 
17 42 




0.05987(45) 
0.05606(47) 
0.05106(49) 
n 0/17^/1 

0.04493(57) 
0.05536(47) 
0.05350(49) 
0.05251(50) 
0.05353(48) 
0.05170(50) 
0.05070(51) 




20 56 
20 56 
20 56 

90 ^« 
ZAj 00 

20 56 
20 56 
20 56 
20 56 
20 56 
20 56 
20 56 


0.04605(67) 
0.04324(72) 
0.03914(80) 

0.03508(89) 
0.04257(73) 
0.04152(76) 
0.04086(80) 
0.04119(76) 
0.04014(78) 
0.03945(81) 


TABLE XXXI. Pseudoscalar meson decay constant Z K fpg. 






/3=5.90 




,9=6.10 


P=6.25 


0=6.47 


u 

fK 


(mx-input) 
(m</>-input) 


0.06882(92) 
0.08343(61) 
0.08781(74) 


0.05201(90) 
0.06216(59) 
0.06486(62) 


0.04172(67) 
0.04970(55) 
0.05170(60) 


0.03304(101) 
0.03913(85) 
0.04074(90) 




TABLE XXXII. Pseudoscalar meson decay constants 


in units of GeV. 








/3=5.90 




/3=6.10 


P=6.25 


/3=6.47 


fK 


(mK-input) 
(ra^-input) 


0.1220(20) 
0.1479(14) 
0.1557(11) 




0.1227(24) 
0.1466(16) 
0.1530(14) 


0.1198(23) 
0.1427(19) 
0.1485(17) 


0.1233(44) 
0.1461(38) 
0.1521(34) 
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TABLE XXXIII. Renormalization constant Zy determined by a non-perturbative method. Ranges of the fit shown at the 
bottom are common to all combinations of quark masses. 





a c on 

p — o.yu 


R a 1 n 

p — D. 1U 


r c\ op: 

p — o.zo 


R a A7 


SjSj 


n c;qp;qq^ oci\ 
vj.ooOovy Zaj) 


u.oyyuy^ zuj 


n fiQQ/i 1 z' oq\ 

U.0OO41^ Zo) 


u.o / / y^ zo j 


S2S2 


n K1 A ( OA \ 


u.oyoo i \ n ) 


U.OoozO^ Ol j 


n ^7/inQ/' Qn^ 

U.O / 4Uo^ OU j 




U.OoooD^ o4 ) 


U.Oy/Dl^ 4Z J 


U.OOZOO^ 04J 


n ^7/infi/' /1Q^ 

U.D / 4UD^ 4y J 


U2U2 


n ^0A(\{ ^0~\ 
\J.ooZ l ^±yjy Oz J 


u.oyoyz^ /uj 


U.OOUUO^lUo j 


U.D / oD 1 ^ 1 O j 


U3U3 


0.53098( 77) 


0.59415(111) 


0.62637(194) 


0.67272(126) 


SlMl 


0.53480( 25) 


0.59838( 28) 


0.63303( 33) 


0.67384( 33) 


S\U2 


0.53400( 29) 


0.59771( 34) 


0.63244( 43) 


0.67353( 39) 


S1U3 


0.53340( 33) 


0.59720( 41) 


0.63173( 54) 


0.67294( 48) 


S2W1 


0.53447( 28) 


0.59816( 33) 


0.63288( 40) 


0.67400( 38) 


S2U2 


0.53368( 33) 


0.59744( 41) 


0.63218( 52) 


0.67371( 45) 


S2U3 


0.53309( 38) 


0.59688( 49) 


0.63134( 66) 


0.67312( 55) 


range 


9 17 


14 22 


20 35 


20 40 




TABLE XXXIV. 


Renormalized vector 


meson decay constant Fya. 






p — o.yu 


/Q — ft in 
p — 0. 1U 


/Q — ft oc: 

p — o.zo 


R a A7 


SlSl 


u.izoyu^ 00 j 


fl flQ/l 7fit A 


ft f\7RA 7t KO\ 


U.UOo4o^ D / ) 


S2S2 


U.lzzlll Do ) 


A fl0O77/ KK\ 


n nv/ioft/ 
U.U/4yo( 00 ) 


U.UODZD(^ OOj 


U\U\ 


Oil QS2V/' 7f<\ 

u. 1 iyo / ^ iO) 


u.uyuo / ^ Dy j 




n CiKfinfit S7^ 

U.UODUD^ / ) 


U2U2 


0.11864( 79) 


0.08960( 69) 


0.07211( 79) 


0.05571(120) 


U3U3 


0.11795( 99) 


0.08912( 88) 


0.07183( 90) 


0.05556(118) 


SlWl 


0.12147( 62) 


0.09233( 55) 


0.07441( 61) 


0.05695( 73) 


S\U 2 


0.12054( 69) 


0.09136( 62) 


0.07397( 56) 


0.05664( 81) 


S\U3 


0.12025( 75) 


0.09106( 69) 


0.07362( 63) 


0.05675( 96) 


S2U1 


0.12081( 67) 


0.09151( 60) 


0.07380( 66) 


0.05662( 78) 


S2U2 


0.12005( 75) 


0.09061( 69) 


0.07345( 61) 


0.05642( 88) 


S2U3 


0.11978( 83) 


0.09029( 77) 


0.07317( 68) 


0.05651(105) 


TABLE XXXV. Vector meson decay constant Fv at physical points in units of GeV. 




P = 5.90 


P = 6.10 


P = 6.25 


P = 6.47 


F P 


0.2256(20) 


0.2227(21) 


0.2166(29) 


0.2163(60) 




0.2416(19) 


0.2430(19) 


0.2360(24) 


0.2321(47) 
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TABLE XXXVI. Pseudoscalar meson masses interpolated or slightly extrapolated to mps/my = 0.75(si), 0.7(s2), 
0.6(mi), 0.5(«2) and 0.4(m3) and normalized by m^?' 75 'a. Parameters and X 1 /dj of their linear continuum extrapolations 
ma — moa + c ■ m^' 75 ^a are also given. 



K1K2 


= 5.90 


= 6.10 


= 6.25 


= 6.47 


moa 


c 


X 2 /df 


Si Si 


0.7500( 0) 


0.7500( 0) 


0.7500( 0) 


0.7500( 0) 


0.7500( 0) 


0.000( 0) 


0.00 


S2S2 


0.6635( 5) 


0.6623( 5) 


0.6624( 6) 


0.6606(16) 


0.6595(15) 


0.008( 4) 


0.38 




0.5282( 8) 


0.5266( 9) 


0.5265(12) 


0.5253(25) 


0.5229(27) 


0.010( 6) 


0.09 


U2U2 


0.4192(11) 


0.4179(12) 


0.4176(16) 


0.4179(31) 


0.4153(36) 


0.007( 9) 


0.07 


W3W3 


0.3243(13) 


0.3235(14) 


0.3231(18) 


0.3247(35) 


0.3223(41) 


0.004(10) 


0.14 


S1W1 


0.6482( 6) 


0.6460( 6) 


0.6462( 8) 


0.6439(29) 


0.6415(20) 


0.013( 5) 


0.66 


S1M2 


0.6074( 7) 


0.6047( 7) 


0.6046( 9) 


0.6029(31) 


0.5989(22) 


0.016( 5) 


0.51 


S1U3 


0.5790( 7) 


0.5757( 7) 


0.5752( 9) 


0.5740(33) 


0.5677(24) 


0.022( 6) 


0.43 


S2UI 


0.5998( 8) 


0.5971( 8) 


0.5974(10) 


0.5941(34) 


0.5915(26) 


0.016( 6) 


0.65 


S2U2 


0.5556( 8) 


0.5525( 8) 


0.5527(11) 


0.5500(35) 


0.5460(28) 


0.018( 7) 


0.58 


S2U3 


0.5244( 9) 


0.5207( 8) 


0.5206(11) 


0.5185(36) 


0.5124(28) 


0.023( 7) 


0.59 



TABLE XXXVII. The same as table XXXVI for vector mesons. Values of o at each are also shown 



K1K2 


= 5.90 


= 6.10 


= 6.25 


= 6.47 


moa 


c 


x 2 /df 


SlSl 


1.0000( 0) 


1.0000( 0) 


1.0000( 0) 


1.0000( 0) 


1.0000( 0) 


0.000( 0) 


0.00 


S2S2 


0.9478( 6) 


0.9462( 7) 


0.9463( 9) 


0.9437(22) 


0.9422(22) 


0.011( 5) 


0.38 


MlWl 


0.8803(14) 


0.8777(15) 


0.8776(20) 


0.8755(42) 


0.8716(46) 


0.017(11) 


0.09 


U2U2 


0.8383(23) 


0.8358(24) 


0.8353(31) 


0.8358(62) 


0.8306(72) 


0.015(17) 


0.07 


U3U3 


0.8108(33) 


0.8087(34) 


0.8078(45) 


0.8118(87) 


0.8056(103) 


0.009(25) 


0.14 


S1W1 


0.9400( 7) 


0.9386( 7) 


0.9386( 9) 


0.9372(21) 


0.9354(22) 


0.009( 5) 


0.13 


S1M2 


0.9187(10) 


0.9173(10) 


0.9171(14) 


0.9167(28) 


0.9142(32) 


0.009( 7) 


0.05 


S1U3 


0.9046(14) 


0.9033(14) 


0.9030(19) 


0.9037(36) 


0.9008(42) 


0.007(10) 


0.08 


S2M1 


0.9140(10) 


0.9118(11) 


0.9118(14) 


0.9094(31) 


0.9067(32) 


0.014( 8) 


0.17 


S2U2 


0.8928(13) 


0.8906(14) 


0.8905(18) 


0.8891(38) 


0.8856(42) 


0.014(10) 


0.07 


S2W3 


0.8788(16) 


0.8767(17) 


0.8764(23) 


0.8763(45) 


0.8723(52) 


0.012(12) 


0.06 


(0.V5) 

rriy a 


0.5086(12) 


0.3874(11) 


0.3212(10) 


0.2447(14) 
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TABLE XXXVIII. The same as table 



XXXVI 



for E-like octet baryons. 



K1K2 


= 5.90 


= 6.10 


= 6.25 


= 6.47 


mo a 


c 


X 2 /df 


SlSlSl 


1.5385(23) 


1.5175(23) 


1.5076(26) 


1.4983(45) 


1.4568(62) 


0.159(15) 


0.31 


S2S2S2 


1.4424(24) 


1.4160(26) 


1.4032(30) 


1.3938(55) 


1.3398(70) 


0.200(17) 


0.66 


MlWlMl 


1.2988(30) 


1.2672(33) 


1.2517(39) 


1.2471(78) 


1.1788(92) 


0.234(22) 


1.44 


U2U2U2 


1.1955(39) 


1.1618(43) 


1.1452(48) 


1.1458(102) 


1.0701(117) 


0.243(28) 


1.73 


U3U3U3 


1.1195(49) 


1.0848(53) 


1.0679(56) 


1.0722(123) 


0.9928(141) 


0.245(34) 


1.71 


MlMlSl 


1.3970(26) 


1.3704(29) 


1.3583(33) 


1.3521(62) 


1.2973(77) 


0.194(18) 


1.08 


U2U2S1 


1.3430(28) 


1.3151(33) 


1.3023(37) 


1.2991(71) 


1.2400(86) 


0.200(21) 


1.53 


U3U3S1 


1.3067(31) 


1.2782(37) 


1.2648(40) 


1.2638(77) 


1.2023(94) 


0.202(22) 


1.75 


SlSlMl 


1.4562(24) 


1.4320(26) 


1.4168(29) 


1.4112(52) 


1.3574(68) 


0.193(16) 


1.21 


S1S1U2 


1.4221(24) 


1.3971(27) 


1.3806(30) 


1.3766(55) 


1.3196(71) 


0.200(17) 


1.70 


S1S1U3 


1.3969(25) 


1.3717(27) 


1.3545(31) 


1.3515(57) 


1.2931(73) 


0.203(18) 


1.97 


H1U1S2 


1.3564(27) 


1.3276(31) 


1.3137(35) 


1.3080(68) 


1.2471(83) 


0.213(20) 


1.22 


M2M2S2 


1.2999(31) 


1.2697(35) 


1.2552(39) 


1.2524(78) 


1.1874(93) 


0.219(22) 


1.60 


M3M3S2 


1.2618(34) 


1.2308(39) 


1.2160(43) 


1.2153(85) 


1.1477(102) 


0.221(24) 


1.78 


S2S2U1 


1.3900(25) 


1.3622(28) 


1.3469(33) 


1.3404(62) 


1.2804(76) 


0.214(18) 


1.13 


S2S2U2 


1.3533(27) 


1.3247(29) 


1.3081(34) 


1.3033(65) 


1.2399(80) 


0.221(19) 


1.46 


S2S2U3 


1.3263(28) 


1.2975(30) 


1.2803(35) 


1.2764(68) 


1.2115(83) 


0.224(20) 


1.62 



TABLE XXXIX. The same as table XXXVI for A-like octet baryons 



Kl«2 


= 5.90 


= 6.10 


= 6.25 


= 6.47 


moa 


c 


x 2 /df 


MlUlSi 


1.3732(26) 


1.3451(28) 


1.3286(33) 


1.3229(62) 


1.2609(77) 


0.219(18) 


1.31 


U2U2S1 


1.3014(29) 


1.2718(31) 


1.2541(36) 


1.2509(72) 


1.1836(87) 


0.230(21) 


1.68 


U3U3S1 


1.2476(33) 


1.2174(34) 


1.1993(39) 


1.1978(80) 


1.1281(96) 


0.233(23) 


1.75 


S1S1U1 


1.4737(24) 


1.4504(25) 


1.4408(29) 


1.4313(53) 


1.3859(68) 


0.171(17) 


0.53 


S1S1U2 


1.4510(25) 


1.4270(27) 


1.4175(30) 


1.4085(56) 


1.3618(71) 


0.174(17) 


0.64 


S1S1U3 


1.4367(25) 


1.4123(29) 


1.4028(31) 


1.3941(59) 


1.3467(73) 


0.175(18) 


0.66 


U1U1S2 


1.3410(27) 


1.3113(30) 


1.2953(35) 


1.2897(68) 


1.2250(83) 


0.226(20) 


1.30 


U2U2S2 


1.2677(31) 


1.2365(34) 


1.2192(39) 


1.2164(80) 


1.1463(94) 


0.236(23) 


1.63 


U3U3S2 


1.2129(36) 


1.1811(38) 


1.1634(43) 


1.1624(90) 


1.0898(105) 


0.239(25) 


1.66 


S2S2U1 


1.4032(25) 


1.3759(28) 


1.3632(32) 


1.3550(61) 


1.2990(75) 


0.203(18) 


0.86 


S2S2U2 


1.3790(26) 


1.3509(30) 


1.3384(34) 


1.3309(64) 


1.2735(79) 


0.205(19) 


1.01 


S2S2U3 


1.3636(27) 


1.3352(31) 


1.3227(35) 


1.3156(66) 


1.2574(81) 


0.207(19) 


1.07 
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TABLE XL. The same as table 



XXXVI 



for decuplet baryons. 



K1K2 


P = 5.90 


P = 6.10 


P = 6.25 


= 6.47 


mo a 


c 


X 2 /df 


SlSlSl 


1.6918(29) 


1.6712(29) 


1.6639(35) 


1.6718(77) 


1.6244(87) 


0.129(21) 


3.02 


S2S2S2 


1.6164(33) 


1.5933(34) 


1.5851(41) 


1.5931(83) 


1.5415(98) 


0.143(23) 


2.94 




1.5205(40) 


1.4948(43) 


1.4845(50) 


1.4969(100) 


1.4374(120) 


0.158(29) 


3.05 


U2U2U2 


1.4629(48) 


1.4354(50) 


1.4227(59) 


1.4402(122) 


1.3725(143) 


0.172(34) 


3.04 


U3U3U3 


1.4271(58) 


1.3979(60) 


1.3827(69) 


1.4051(151) 


1.3286(172) 


0.188(41) 


2.75 


MlMlSl 


1.5768(36) 


1.5529(37) 


1.5435(44) 


1.5553(89) 


1.4999(106) 


0.146(25) 


3.45 


U2U2S1 


1.5372(40) 


1.5123(42) 


1.5012(48) 


1.5176(99) 


1.4569(118) 


0.153(28) 


3.88 


U3U3S1 


1.5118(45) 


1.4859(48) 


1.4730(53) 


1.4943(111) 


1.4278(131) 


0.159(32) 


4.22 


SlSlMl 


1.6339(32) 


1.6117(33) 


1.6033(39) 


1.6136(81) 


1.5622(94) 


0.136(23) 


3.53 


S1S1U2 


1.6135(33) 


1.5909(35) 


1.5816(40) 


1.5948(83) 


1.5408(98) 


0.138(24) 


4.08 


S1S1U3 


1.6000(35) 


1.5770(37) 


1.5667(42) 


1.5832(87) 


1.5265(103) 


0.139(25) 


4.69 


H1U1S2 


1.5522(37) 


1.5274(40) 


1.5177(47) 


1.5290(93) 


1.4722(111) 


0.152(27) 


3.14 


M2M2S2 


1.5129(42) 


1.4870(44) 


1.4757(51) 


1.4913(104) 


1.4291(124) 


0.159(30) 


3.46 


M3M3S2 


1.4877(47) 


1.4609(50) 


1.4478(56) 


1.4680(117) 


1.4001(138) 


0.167(33) 


3.71 


S2S2U1 


1.5841(35) 


1.5602(37) 


1.5513(44) 


1.5611(88) 


1.5069(104) 


0.147(25) 


3.11 


S2S2U2 


1.5641(37) 


1.5396(39) 


1.5299(45) 


1.5422(91) 


1.4855(109) 


0.150(26) 


3.46 


S2S2U3 


1.5508(39) 


1.5260(41) 


1.5153(47) 


1.5306(95) 


1.4711(114) 


0.152(27) 


3.88 



TABLE XLI. Quark masses m q a in the MS scheme at fj, = 2 GeV interpolated or slightly extrapolated to 
mps/niv = 0.75(si), 0.7(s2), 0.6(iti), 0.5(112) and 0. 4(113) and normalized by m^?' 75 'a together with parameters and xV^f °f 
their linear continuum extrapolations m q a — m q fia + c- my' 7 ^ a. Lattice spacings at each /3 determined from the physical value 
of my 7 '' are also given. 



K1K2 


P = 5.90 


P = 6.10 


P = 6.25 


P = 6.47 


m 9i oa 


c 


x 2 /df 


Si 


0.1164( 3) 


0.1185( 4) 


0.1203( 5) 


0.1216(10) 


0.1265(11) 


-0.020( 3) 


0.28 


S2 


0.0910( 3) 


0.0927( 4) 


0.0943( 5) 


0.0949( 9) 


0.0992(10) 


-0.016( 2) 


0.51 


Ul 


0.0568( 3) 


0.0581( 3) 


0.0593( 4) 


0.0598( 8) 


0.0630( 9) 


-0.012( 2) 


0.40 


U2 


0.0347( 3) 


0.0357( 3) 


0.0365( 4) 


0.0370( 8) 


0.0393( 9) 


-0.009( 2) 


0.15 


U3 


0.0199( 2) 


0.0205( 2) 


0.0209( 4) 


0.0215( 7) 


0.0226( 8) 


-0.005( 2) 


0.05 


a-' (GeV) 


1.929(36) 


2.532(47) 


3.054(57) 


4.009(78) 
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TABLE XLII. Parameters of QxPT fits made in the continuum limit. Results are given in units of m v ' of the continuum 
theory and in physical units (GeV,GeV _1 ). We also give parameters obtained from an extrapolation of the parameters at finite 
f3. 





X 2 /N df 


Par am. 


in rriy' 7 ^ unit 


in Phys. unit 


extrap.(in Phys. unit) 


PS 


0.2 


A 
B 
5 


f.72(23) 
2.17(91) 
0.16(10) 


1.69(23) 


2.43(4) 
1.69(9) 


Vector 


0.0 


m° v 


0.791(29) 


0.776(11) 


0.777(10) 






Ci 


0.548(79) 


0.559(93) 


0.560(72) 






Cl/2 


-0.133(98) 




-0.136(82) 




0.0 


IliQ 


8771421 


861 (20) 


863^24^ 






b F 


-0.248(66) 


-0.253(67) 


-0.244(36) 






bo 


0.040(44) 


0.041(45) 


0.040(18) 






F 


0.290(71) 




0.280(36) 






D 


0.434(69) 




0.416(35) 


decuplet 


0.0 


m° D 


1.265(72) 


1.241(74) 


1.251(43) 






c 


0.220(77) 


0.224(78) 


0.234(28) 






H 2 


0.06(99) 




0.16(40) 






C 2 


0.3(3.0) 




-0.38(77) 



TABLE XL11L Hadron spectrum from Q^PT chiral fits in the continuum limit. The scale m v of the continuum theory 
is also given. See text for details. Deviation and its statistical significance are relative to our final result (data in fVl)- 





vtik input 




input 






mass(GeV) 


deviation 


mass(GeV) 


deviation 




0.981(18) 


0.981(18) 


K 






0.559( 6) 


1.0% 0.6a 


K* 


0.856( 9) 


-0.3% 0.3a 


0.889( 4) 


0.0% 0.0a 




0.951(12) 


-0.6% 0.4a 






N 


0.876(23) 


-0.2% 0.1a 


0.876(23) 


-0.2% 0.1a 


A 


1.011(17) 


-0.8% 0.4a 


1.057(12) 


-0.3% 0.3a 


E 


1.109(16) 


-0.7% 0.4a 


1.173(10) 


-0.2% 0.2a 




1.189(13) 


-1.0% 0.7a 


1.285( 8) 


-0.3% 0.5a 


A 


1.251(48) 


-0.5% 0.2a 


1.251(48) 


-0.5% 0.2a 


E* 


1.353(32) 


-0.4% 0.2a 


1.385(31) 


-0.2% 0.1a 




1.453(27) 


-0.4% 0.2a 


1.517(25) 


0.0% 0.0a 


n 


1.551(20) 


-0.6% 0.4a 


1.647(14) 


0.0% 0.0a 
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TABLE XLIV. 


Smearing parameters 


for quark matrix inversion. Parameters A and B 


appear in the smearing function 


S(n) = Acxp(-B 


n |). The number of iterations necessary for point (P) and smeared (S) sources are also given. 








P = 5.90 (rp = 10~ 12 


, r% = 10 7 ) 






K 


0.15660 


0.15740 


0.15830 


0.15890 


u. loyzu 


4 


1.1427 


1.1408 


1.1387 


1.1373 


1 1 Odd 
l.lODD 


JD 
JD 


0.3290 


0.3268 


0.3242 


0.3226 


A QOI <7 

U.ozl ( 


#lter(FJ 


220±8 


307±13 


495±26 


70H 1 1 ta 

780±170 


1 1 on_Loon 


#lter(b) 


201±6 


279±10 


A a O A oo 


750±170 


i nori-Lonn 






P = 6.10 (rp = 10~ LZ 


2 -i a-7\ 

,rs = W ') 






K 


0.15280 


0.15340 


0.15400 


0.15440 


O 1 ^A«f\ 

u. id^ou 


A 
f\ 


1.0884 


1.0863 


1.0843 


1.0830 


l.Uozo 


JD 
JD 


0.2361 


0.2339 


0.2317 


0.2303 


a oont? 
U-zzyD 


#lter(FJ 


259±6 


373±9 


644±21 


1170±60 


1 *7 4f\ 1 OOfl 

1 / 4U±zoU 


#lter(bj 


261±4 


378±7 


661±18 


1220±60 


i Q/i a_i_o e: a 






P = 6.25 (rp = 10~ 12 


7 1 n — 7 \ 

, = 10 ') 








0.15075 


0.15115 


0.15165 


0.15200 


n i coon 
U. 10ZZU 


4 
71 


1.06304 


1.06151 


1.05961 


1.05829 


1 A 

l.UO/04 


7~> 
£> 


0.19336 


0.19175 


0.18974 


0.18835 


A 1 07C C 


#lter(F) 


275±6 


359±8 


586±15 


1080±40 


1 AAA_I_ 1 OA 

iyuu±izu 


#lter(b) 


298±3 


387±5 


633±11 


1170±40 


O 1 Q A_l_ 1 OA 






P = 6.47 (rp = 10" 14 


.rS = 10-') 






K 


0.14855 


0.14885 


0.14925 


0.14945 


0.14960 


A 


1.02817 


1.02688 


1.02517 


1.02431 


1.02368 


B 


0.13484 


0.13351 


0.13174 


0.13085 


0.13019 


#Iter(P) 


427±8 


561±11 


987±63 


1600±50 


3060±160 


#Iter(S) 


405±3 


532±6 


940±60 


1520±40 


2930±140 



TABLE XLV. Values of Xfull f° r degenerate (DG) and non-degenerate (ND) pseudoscalar meson mass fits. 







DG 




ND (ki = si) 




ND («i = si) 


p 


N df 


Xfull 


N d f 


2 

Xfull 


N df 


2 

Xfull 


5.90 


85 


79(21) 


68 


61(17) 


68 


61(18) 


6.10 


100 


119(26) 


80 


92(22) 


80 


99(23) 


6.25 


40 


47(15) 


32 


29(12) 


32 


33(13) 


6.47 


70 


172(66) 


56 


85(36) 


56 


86(35) 
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TABLE XLVI. Chiral fit functions and Xext/Ndf for pseudoscalar meson masses. Fit 1 - Fit 5 are for degenerate cases, while 



Fit In - 
rows. 


Fit 5n are for non-degenerate cases with 


m s = (1/si - 


1/k c )/2 for upper rows 


and m s — 


(l/s2 — l/« c )/2 for lower 




fit function 


parameters 


N df 


P = 5.90 


Xext 

13 = 6.10 


/N df 

f3 = 6.25 


13 = 6.47 


Fit 1 
Fit 2 
Fit 3 
Fit 4 

Fit 5 


'2 O A 

m PS = 2 Am 
m 2 PS = 2Am + 45m 2 
m 2 PS = 2Am + 45m 2 + 8Cm 3 
mis = 2Am\\ - 5(\n{f A m/A) + 1)] 
mp s = 2Am[l - 8(ln(f A m/A) + 1)] 
+45m 2 


k c ,A,B 
k c ,A,B,C 
He, A, 8 

k c ,A,S,B 


3 
2 
1 
2 

1 


1290(100) 
417( 43) 
182( 43) 

1018( 85) 

116( 31) 


1117( 92) 
172( 28) 
125( 54) 
606( 61) 

71 ( 35) 


942( 82) 
67( 28) 
100( 55) 
376( 55) 

57( 31) 


630(170) 
59( 30) 
83( 48) 

285( 92) 

26( 26) 


Fit In 
Fit 2n 

Fit 5n 


m 2 PS = A a m s + Am 

m 2 PS = A(m s + m) + B(m s + m) 2 

+C(m s - m) 2 
m 2 PS = A(m s + m)[l — 8{\n{f A m a /A) 

+m/(m — m s ) ln(m/m s ))] 

+B(m s +m) 2 


A s , A 
A,B,C 

A,S,B 


2 
1 

1 


630(110) 
381( 84) 
91( 26) 
128( 35) 
147( 40) 

80( 50) 


455(55) 
314(45) 
20( 12) 
36(17) 
46(19) 

62( 24) 


430( 50) 
339( 41) 
2(3) 
4( 4) 
12(9) 

14( 10) 


246( 67) 
233( 64) 
0( 1) 
1(4) 
4( 9) 

8( 13) 



TABLE XLVII. Bare AWI quark masses 2mf WI ' local< - ^ obtained at the simulation points. Here, 2m q stands for m qf +rn, 



The fit range [t m in,t m ax] is also given. 





f3 = 5.90 


P = 6.10 


P = 6.25 


P = 6.47 


K1K2 


range 


value 


range 


value 


range 


value 


range 


value 


Si Si 


5 25 


0.091309( 80) 


8 32 


0.072038( 59) 


7 39 


0.063377(74) 


10 53 


0.048770( 80) 


S2S 2 


5 25 


0. 067984 ( 74) 


8 32 


0.051800( 53) 


7 39 


0.048763(67) 


10 53 


0.036887( 72) 


MlMl 


5 25 


0.042573( 68) 


8 32 


0.031987( 50) 


7 39 


0.030771(57) 


10 53 


0.021232( 60) 


M2M2 


5 25 


0.025941( 69) 


8 32 


0.018921( 56) 


7 39 


0.018359(52) 


10 53 


0.013501( 58) 


«3"3 


5 25 


0.017144(142) 


8 32 


0.012114(100) 


7 39 


0.011225(60) 


10 53 


0.007701( 62) 


SlWl 


5 25 


0.066713( 75) 


8 32 


0.051875( 54) 


7 39 


0.046989(66) 


10 53 


0.034941( 73) 


Sl«2 


5 25 


0.058270( 75) 


8 32 


0.045274( 54) 


7 39 


0.040706(64) 


10 53 


0.031035( 73) 


S1U3 


5 25 


0.054006( 78) 


8 32 


0.041953( 57) 


7 39 


0.037117(65) 


10 53 


0.028115( 78) 


S 2 Wl 


5 25 


0.055217( 71) 


8 32 


0.041860( 51) 


7 39 


0.039740(62) 


10 53 


0.029039( 67) 


S2U 2 


5 25 


0.046842( 71) 


8 32 


0.035299( 52) 


7 39 


0.033486(59) 


10 53 


0.025147( 66) 


S2U3 


5 25 


0.042594( 75) 


8 32 


0.031982( 55) 


7 39 


0.029911(60) 


10 53 


0.022235( 70) 
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FIG. 1. Chiral extrapolations for (a) pseudoscalar meson, (b) vector meson, (c) octet baryon and (d) decuplet baryon at 
f3 = 5.9. The QxPT and polynomial chiral fits are shown by solid and dashed lines. The insets are expanded dis play s for 



degenerate cases. In panel (a) AWI quark mass, mf WI and linear chiral extrapolation are given (discussed in Sec. VII). In 



panels (c) and (d), we give data only for combinations of (s\,Ui,Ui) and (t 
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FIG. 5. Continuum extrapolation of light hadron masses 
from mx-input. Filled symbols and solid lines are the results 
from the QxPT chiral fits, while open symbols and dashed 
lines are from the polynomial fits. Experimental values are 
shown by the stars. 




FIG. 6. Quenched light hadron spectrum compared with 
experiment. The statistical error and the sum of the statisti- 
cal and systematic errors are indicated. 
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FIG. 7. Hyperfine splitting of mesons normalized by tuk- 
Diamonds represent the experimental points corresponding to 
(mps,mv) = (mn,m,p) and (m,K,m* K ), where the former is 
the input. 
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FIG. 8. J parameter. The star represents the experimental 
value. 
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FIG. 9. Comparison of the spectra from method A (black 
circles) and method B (open circles), ma is taken as input. 
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FIG. 10. Edinburgh plot in the continuum limit. The stars 
represent experimental values. Dashed curves illustrate the 
phenomenological mass formulae by Ono po| . 
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FIG. 11. The averaged mass of the up and down quarks 
and its continuum extrapolation. The leftmost point is the 
value extrapolated to the continuum limit. 
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FIG. 12. The strange quark masses and their continuum 
extrapolation. The leftmost points are the values extrapo- 
lated to the continuum limit. 
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FIG. 13. Continuum extrapolations of pseudoscalar meson 
decay constants / w and fx (jtik- input). Large squares at 
a — represent experimental values. 
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FIG. 14. Continuum extrapolation of /a'//tt — 1 from the 
mic-input. The square at a = represents the experimental 
value. 
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FIG. 15. F p (top panel) and (bottom panel) as a func- 
tion of a. Fits for the continuum extrapolation are also shown. 
Solid symbols are for non-perturbatively renormalized decay 
constants, and open symbols are for decay constants renor- 
malized by tadpole-improved one-loop perturbation theory. 
Stars represent experimental values. 
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FIG. 16. Typical effective masses obtained with various 
combinations ol quark sources. 
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FIG. 17. Effective mass plots of degenerate octet baryons 
at mps/mv « 0.75 (top), mps/mv ~ 0.6 (middle), and 
mps/mv ~ 0.4 (bottom) at (3 = 5.90. 
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FIG. 18. The same as Fig.fj]] at /3 = 6.10. 
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FIG. 19. The same as Fig.0 at f3 = 6.25. 
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FIG. 20. The same as Fig.fj]] at (3 = 6.47. 
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FIG. 22. The same as Fig. El at (3 = 6.47. 
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FIG. 21. Effective mass plots for twice the AWI quark 
mass, 2m^ wi ^ ', at (3 = 5.90 for the degenerate cases, cor- 
responding to mps/mv ~ 0.75 (top), 0.6 (middle) and 0.4 
(bottom). 
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FIG. 23. Deviations from 

the chiral relation y = rn 2 PSl2 /^ rn i WI + rn,2 WI ) — constant. 
The horizontal axis is x = mf WI + m^ wl . Filled and open 
symbols represent degenerate and non-degenerate quark mass 
cases, respectively. The error from the lattice spacing is not 
included. 
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FIG. 24. Test of quenched chiral logarithms for pseu- 
doscalar meson masses. 
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FIG. 25. Ratio test proposed in Ref. [§. Data at /3 = 5.90 
are connected to the point (0.0, 1.0) to guide the eyes. 
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FIG. 26. Ratio test of quenched chiral logarithms for pseu- 
doscalar meson decay constants. 
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FIG. 27. Degenerate hadron masses vs. mps at /3=5.9. 
The leftmost points are values extrapolated to the chiral limit, 
and the second ones from the left are those at the physical 
point. Fits from two types of chiral extrapolations based on 
QxPT are shown. See text for details. 
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FIG. 28. Light hadron spectrum from polynomial chiral fits. 
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FIG. 29. Comparison of our hadron masses (filled sym- 
bols) with those of Ref. || (open symbols). Open triangles 
(squares) are data obtained on lattices with L s a w 2.3 fm 
(L s a « 3.4 fm). Open symbols at a « represent data on fi- 
nite lattice (right) and those after finite size corrections (left). 
Experimental values are shown by stars. 
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FIG. 30. The ratio y = {mX WI +mT 
as a function of x = ra^ WI + m<t WI . Filled and open sym- 
bols are for degenerate and non-degenerate cases, respectively. 
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FIG. 31. Chiral extrapolations of the quark mass based on 
the axial- vector Ward identity at /3 = 5.90. 
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FIG. 32. Comparison of the light quark masses determined 
from various chiral fits. The VWI (AWI) quark masses are 
shown by filled (open) symbols. Circles are results obtained 
from our main analysis, while triangles are from the alter- 
native chiral fits discussed in the text. Open diamonds are 
obtained from fits to m PS as a function of m^ WI . 
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FIG. 35. Comparison of the u, d quark mass. Circles are 
the VWI quark masses. Triangles (squares) are the AWI 
quark masses derived from the extended (local) axial vector 
current. Filled (open) symbols indicate that masses are cal- 
culated using ap (*Mg-)- The leftmost data shows the result 
from a combined fit with the statistical error and the sum of 



statistical and systematic errors. See Sec. VII for detail 
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FIG. 38. A^g vs a. 



0.0 



0.2 



0.4 



a (GeV 1 ) 



FIG. 39. yfa/rrip and l/rom p 



0.6 









, -Cr ' 












• VCT/m p 




1/(r,m,) 



0.6 



0.49 
0.48 
0.47 
0.46 
0.45 
0.44 
0.43 
0.43 
0.42 
0.41 
0.40 
0.39 
0.38 
0.41 
0.40 
0.39 
0.38 
0.37 
0.36 



fflS 



10 15 20 25 30 35 
t 



FIG. 40. Examples of fits for Eq. @ at j3 = 6.10 in units 
of 10 . Three panels show the data for mps/mv ~ 0.75 (top), 
0.6 (middle) and 0.4 (bottom). 
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FIG. 41. Pseudoscalar meson decay constants versus \/k 
at j3 = 5.90. Solid curves show chiral fits with quadratic 
polynomial (linear) for the degenerate (non-degenerate) case. 
Dashed line is for a linear fit to the degenerate case. 
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FIG. 42. Ratios (V c V L )/(V L V L ) and fits to obtain Z v at 
(3 — 6.25. Parameters are mps/mv ~ 0.75 (top), mps/my « 
0.6 (middle), and mps/mv ~ 0.4 (bottom) at (5 = 5.90. 
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FIG. 43. Fya vs. 1/k at /3 = 5.90 with chiral extrapola- 
tions. Non-degenerate data are slightly shifted in x for clarity. 
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FIG. 44. Degenerate masses and chiral fits in the contin- 
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FIG. 45. Octet baryon mass in the chiral limit versus the 
lattice spacing. The leftmost point represents the value in 
the continuum limit but determined from the chiral fit in the 
continuum limit. 
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FIG. 46. Ratio test for Q^PT formulae for degenerate vec- 
tor mesons. 
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FIG. 48. The same as Fig. |f6j, but for degenerate decuplet 
baryons. 
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FIG. 47. The same as Fig. [l6| but for non-degenerate vec- 
tor mesons. 
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FIG. 49. The same as Fig. but for degenerate octet 
baryons. 
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